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Abstract 

This is the second of a pair of articles on scattering of glue by glue, in which we give the light-cone gauge 
calculation of the one-loop on-shell helicity conserving scattering amplitudes for gluon-gluon scattering 
(neglecting quark loops). The factors in the gluon propagator are regulated by replacing integrals 
with discretized sums omitting the = terms in each sum. We also employ a novel ultraviolet regulator 
that is convenient for the light-cone worldsheet description of planar Feynman diagrams. The helicity 
conserving scattering amplitudes are divergent in the infra-red. The infrared divergences in the elastic 
one-loop amplitude are shown to cancel, in their contribution to cross sections, against ones in the cross 
section for unseen bremsstrahlung giuons. We include here the explicit calculation of the latter, because 
it assumes an unfamiliar form due to the peculiar way discretization of p"*" regulates infrared divergences. 
In resolving the infrared divergences we employ a covariant definition of jets, which allows a transparent 
demonstration of the Lorentz invariance of our final results. Because we use an explicit cutoff of the 
ultraviolet divergences in exactly 4 space-time dimensions, we must introduce explicit counterterms to 
achieve this final covariant result. These counter-terms are polynomials in the external momenta of 
the precise order dictated by power-counting. We discuss the modifications they entail for the light-cone 
worldsheet action that reproduces the "bare" planar diagrams of the giuonic sector of QCD. The simplest 
way to do this is to interpret the QCD string as moving in six space-time dimensions. 
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1 Introduction 



This is the second of a pair of articles on scattering of glue by glue to one loop in the language of the 

lightconc worldsheet [1-3]. The first article (I) [4] worked out the (finite) one loop amplitudes for helicity 
nonconserving processes. In this article we extend the calculation to the helicity conserving case, for which 
it is necessary to deal with ultraviolet and infrared divergences. 

We refer the reader to the introduction of I for the detailed motivation and background for this work. 
Here we briefly mention the highlights. The goal of the program to give a worldsheet description of the 
sum of planar diagrams is to shed light on Field/String duality from the "field" side at weak coupling. This 
is just the perturbation expansion of the field theory. The mapping of the sum of planar diagrams to a 
worldsheet system [2, 3] in essence allows one to read off the worldsheet dynamics for scalar field theory 
and Yang-Mills theory in this weak coupling limit. A serious limitation of these initial articles, however, is 
that they transcribed the "bare" Fcynman diagrams, without including any of the countcrterms necessary to 
maintain gauge invariance. In the context of dimensional regularization this limitation is innocuous, because 
dimensional regularization automatically includes them. However, we want the worldsheet formalism to 
work in four dimensions, and so we must have in the worldsheet action the flexibility to include countcrterms 
that go beyond the initial input Lagrangian. In [5] the ultraviolet structure of (j)^ theory was analyzed on 
the lightcone and it was shown, to all orders in perturbation theory, that two new counterterms in addition 
to those associated with mass, wave function, and coupling renormalization were necessary and sufficient. 
Happily, a local modification of the "bare" worldsheet action allowed for these new terms. 

The aim of I and the present article is to execute the same program for Yang-Mills theory in light-cone 
gauge. Because the corresponding analysis is considerably more complex, we have limited their scope to one 
loop. In I we focused on one-loop helicity-violating amplitudes for which the on-shell tree diagrams vanished. 
As a consequence the one loop amplitudes are finite in both the ultraviolet and infrared. We could therefore 
confirm that the worldsheet description produces the correct known answers without dealing with coUinear 
and soft gluon emission processes. 

In contrast, the helicity- conserving processes studied in the present article display the full infrared di- 
vergence structure of non-abelian gauge theory. Just as in I, our infrared regulator is discretization of the 

integrals omitting the p"*" = terms. But for the one-loop amplitudes we deal with here, this is essen- 
tially equivalent to simply reserving the p+ integrations to last. This is because all of the artificial = 
divergences actually cancel algebraically if the integrands from all diagrams, with momentum routed appro- 
priately, are combined before the loop integrations arc performed [6]. However, the true infrared divergences, 
for which discretization provides a temporary regulator, are not cancelled in this combination, but rather 
are organized to cancel against divergences due to the absorption and emission of real soft and collincar 
gluons according to the Lec-Nauenberg theorem [7]. In addition, of course, these amplitudes have ultraviolet 
divergences which are taken care of by the renormalization program [8,9]. In this article we give a complete 
lightcone gauge calciilation of the scattering of glue by glue, including the processes involving extra gluons 
that resolve the infrared divergences. We work in four dimensions using the worldsheet friendly ultra-violet 
cutoff employed in [4,5]. We organize the Feynman diagrams of the SU{Nc) Yang-Mills theory according 
to 't Hooft's large Nc expansion [1], and we calculate the one-loop planar diagrams surviving the Nc — > oo 
limit. The 't Hooft limit suppresses diagrams with quark loops, so they are not included here. 

We remark here that our non-traditional methods of regulating and dealing with divergences are not 
sheer perversities on our part, but rather they are guided by the desire to fit these calculations into the 
framework of the lightcone worldsheet formalism. We hope the reader will bear with us on this point, and 
we assure her that, although unconventional, we have been extremely careful with the well-known subtleties 
and pitfalls of this difficult subject [10, 11]. 

To keep this paper reasonably self-contained, we repeat two short sections from I that summarize the 
lightcone Feynman rules (Section 2) and some useful identities (Section 3). A brief Section 4 lists all the 
four gluon trees. Then in Section 5 we discuss bremsstrahlung processes. We use a covariant definition 
of jets that proves to be very helpful in achieving nice compact results for these processes. We deal with 
initial state collinear (mass) divergences as in the original Lee-Nauenberg paper [7] , where they are cancelled 
by including extra near collinear gluons in the initial state. This is in contrast to the standard technique. 
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used in analywing jets experimentally, that absorbs the initial state collinear divergences into the initial state 
parton distribution function [11]. Section 6 briefly summarizes the results for triangle and swordfish diagrams 
obtained in [12]. Section 7 contains the meat of this paper, the calculation of box diagrams. The reduction 
procedure developed in [4] for the helicity violating box diagrams is helpful for this. Section 8 describes 
the results of calculating the remaining quartic triangle and double quartic diagrams, and Section 9 finally 
puts everything together, and explains the cancellation of infrared divergences in cross sections. Section 10 
discusses the problem of giving a local worldsheet description of all the counterterms necessary for Lorentz 
invariance. We find this can be simply done by interpreting the "string" dynamics of the worldsheet as 
occurring in six dimensional spacetime. The two added dimensions are holographically generated on the 
"string" side of Field/String duality and play no role on the "field" side. A final Section 11 wraps up the 
paper with a brief look at issues still to resolve in the future. 

The reader who does not wish to follow the technical details of this work may get a glimpse of the main 
results and their impact on the lightcone worldsheet dynamics by simply reading sections 2, 9, 10 and 11. We 
particularly would like to draw his attention to the nice and compact final results for the infrared finite and 
Lorentz covariant probabilities (more precisely, the unnormalized squared amplitudes for jet-jet scattering) 
Eqs. (121), (122) of Section 9.2. The implications of the necessary counterterms for the worldsheet dynamics, 
which is the subject of Section 10, should also be amusing for this reader. 



2 Feynman Rules for Light-cone gauge Yang-Mills 

Here, we use the notation and conventions in Ref. [13], according to which the values of the non- vanishing 
three transverse gluon vertices are: 




P1P2 P1P2 



(1) 




ptP2 
1 2 

The quartic vertices in this helicity basis are given by 



^{ptp'2-Ptpr) = ^K\ 
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' {Pt+Pt? ^ ^ 
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' {Pt+Ptf {Pt+Pt? J ^ ' 

In these expressions, = (p^ + zp^)/-\/2, = (p^ — ip^)/\/2, and = (p° +p|)/-\/2 are momenta entering 
the diagram on leg fc, and g is proportional to the conventional QCD coupling g^. Note that these are light- 
cone gauge {A- = 0) expressions and include the contributions that arise when the longitudinal field Aj^ is 
eliminated from the formalism.^ We also should point out that we are giving these rules in the context of the 
't Hooft's 1/Nc expansion at fixed Ncgl- Then the planar diagrams of the SU{Nc) theory are correctly given 

^Thcsc vertex rules arc convenient for the mixed r = ix^ , p+ , p representation used in the imaginary 2;+ worldsheet formalism 
[2], in which an i from each vertex has been absorbed in each dx^: idx'^ = dr and the propagator is (2p+)^^e^'^P /^P"*". It 
will sometimes be convenient in this paper to return to full Minkowski momentum space p^,p+,p. Then with the vertex rules 
given here and no i's in the momentum space propagators 1 / (p^ — it) with = — 2p'^p~ , each Minkowski loop momentum 
integral should be accompanied by a —i: —id^qM- With a further Wick rotatation to Euclidean space d^qM = id^QE the i's 
would disappear entirely. 
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if wc take g = \J Nc/2. Non-planar diagrams with this definition of g must be aceompanied by appropriate 
powers of l/N^, depending on the number of "handles" in the diagram. We have not spelled the details out 
here, because our focus will be on the planar diagrams in this article. The results we obtain should therefore 
be compared to the limit Nc oo, fixed g^Nc of those in the literature. In making such comparisons, note 
that our definition of g multiplies conventionally defined n-gluon tree amplitudes by a factor Nc^^~^ Nc 
for n = 4, so for each gluon scattering process we remove this factor before comparing to the literature. 



3 K Identities 

As we have seen, the quantities 

Ktj ^ Ptp^-Ptp^ (5) 

play a central role in the cubic Yang-Mills vertex. In fact, wc shall find that the simplest forms of the various 
helicity amplitudes are achieved by expressing them as functions of the K^s. These simple forms are in fact 
identical to those achieved by Parke and Taylor using a bispinor representation of polarization vectors as the 
now famous Parke- Taylor amplitudes [14]. For us the role of the spinor matrix elements in those formulae 
will be played exclusively by K^j and /ST^^-. 

In order to reduce the expressions for the helicity amplitudes to the Parke- Taylor form, we will need 
a number of identities enjoyed by the JC's. For a general n-gluon amplitude wc can form Kij for each 
pair of gluons (ij), where i,j = 1, . . . ,n distinguish the different gluons. By momentum conservation, it is 
immediate that 

E^i^- = 0- (6) 
j 

Prom the fact that K is an anti-symmetric product we have Bianchi-like identities 

PfK^k+PtKt.+ptKi:, = (7) 
K^,K!>^ + K^^K^^+K^^K^, = (8) 

Finally, the most powerful type of identity follows from a very simple calculation 

] ^3 3 ^3 

which seems like a complicated non-linear relation. However when wc arc considering scattering amplitudes, 
the momenta all satisfy p? = so the right side is zero! This identity plays a central role in showing that 
trees with all but one like helicities vanish. (Trees with all like helicity can't even be drawn.) They are also 
crucial for reducing the complexity of the helicity amplitudes that don't vanish. 



4 Summary of Tree Amplitudes 

In this section we simply list the four point trees amplitudes on and off shell obtained in [4]. There arc no 
tree diagrams for A A AA or VV VV polarizations. The off-shell A A AV four-point tree is given by, omitting 
the coupling factor 2g for each vertex. 



Pt 



PiPtpt 



^32^14, 



+ 



-"^ 43-^^21 
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On shell 



(10) 
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The notation liorc is that A denotes an incoming arrow representing helicity +1, while V denotes an outgoing 
arrow representing hehcity —1. 

The only non-zero four point trees are those with two of each helicity. There are two distinct helicity 
patterns. The amplitude for adjacent helicity A A VV is given by 
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When some legs are off-shell, we use the shorthand notation p* = p^/pf to simplify the writing. 

Atree _ ■^^21 ^43 
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The other distinct helicity arrangement for four gluon scattering is alternating helicity A V AV: 

1 



(13) 
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where the quartic vertex contribution has been split between the last terms in each of the square brackets. 
Notice that the second line on the right side is obtained from the first line with the relabeling substitutions 
1^2^3^4^1 and A ^ V ^ A. Furthermore the first line can be obtained from the first line on the 
right of (11) by interchanging 2 3 and multiplying by the factor —1. Thus by inspection we immediately 
obtain the simplifications 
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Here and in the following we use the shorthand notation Pi,i+i = pi +pi+i, with j = 0, 1, 2, 3 and p4 = po- 



5 Gluon Bremsstrahlung and Jets 

The consistent resolution of infrared divergences in loop corrections to scattering amplitudes involves a 
cancellation against corresponding infrared divergences in the cross section for the emission (or absorption) 
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of an extra gluon, whose momentum is either eolhnear with one of the gluons in the core process or "soft" . 
If the core process is scattering of glue by glue, the associated bremsstrahlung amplitudes are 5 gluon 
amplitudes. 

In the context of the large Nc limit it is necessary to combine coherently only the bremsstrahlung diagrams 
with the same cyclic ordering. So, for example, in the diagrams shown in Fig. 1 at A^c = oo it is only necessary 
to square the sum of the two diagrams on each line and combine the results on different lines incoherently. 
Because Nc = oo suppresses nonplanar diagrams, a gluon line attached between two outgoing gluons (as 




Figure 1: The bremsstrahlung diagrams associated with glue-gluc scattering involving leg 4. At Nc = oo 
the sum of the diagrams on each line may be independently squared to give the leading contribution to the 
cross section. Similar pairs of diagrams involving each of the other legs must also be included. 

with the diagrams on the first line of Fig. 1) must be outgoing when Nc = oo^ . Similarly a gluon line 
attached between two incoming gluons must be incoming. On the other hand, a gluon attached between an 
incoming and an outgoing gluon (as with the diagrams on the second line of Fig. 1) may be either incoming 
or outgoing. 

As is well known infrared and collinear divergences are present only when the bremsstrahlung gluon 

attaches to external legs. For example if the brem gluon is collinear with p4, there is a collinear divergence 
in the phase-space integral of the square of the diagrams where the gluon is emitted from or absorbed by leg 
4. Calling the brem gluon's four-momentum fc, for fixed k'^ the collinear point is fe = k'^p^/p^. Then it is 
convenient to write 

fe = A;+^+fe (16) 

Pi 

and examine the phase space integral for in a neighborhood of zero. Here we assume A;"*" = 0(1) so the 
brem gluon is not soft. In effect, rather than measuring a single gluon, we insist that we measure a "jet" 
of total trausNcrso luonK^utuni P % (/. + +pt)Pi/pt within a resolution A [15]. A simple calculation shows 

"^At first glance the reader might think that the distinction between an incoming or outgoing line has nothing to do with 
planarity. But a well-defined large limit only makes sense for physical quantities that are singlets under the gauge group. 
We specify our planar scattering amplitudes by imagining that the external gluons are all attached to a huge connected Wilson 
loop including a single connected portion at late times and another single connected portion at early times. The large Nc limit 
of such a quantity has the properties we describe in this paragraph. 
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that 

\k+\\pt\ 



= 2|fc||p4|(l-cos(9) (17) 



where the overarrow denotes the 3 spatial components of a four- vector and cos 9 = k ■ p4/\k\\p4\. The left 
side is thus a nice measure of the angular size of a jet, so we define the phase space of a jet of resolution A 
by the restriction 



\k+\\pt\ 



~2 



This translates to fe < |fc+|A2/|p+|. 

The amplitudes for the emission of a hard brcm gluon from the right of leg 4 (as in the first diagram on 
the first line of Fig. 1) are given, for the two polarizations, by 



^Brem = ''^9^^ ' (k + pi)^ ' Outgoing helicity (19) 

"^^-^ = -^ h^k^+pt) (fcT^I? ' (20) 

When the brem gluon (with momentum k) is emitted from the left of leg 4, the amplitudes are the same 
except that -K"4,fe appears instead of -K'fe,4. Thus the amplitudes for emission from left and right have opposite 
signs. The amplitudes do not cancel, however, because they have different gauge group structure. As already 
mentioned at Nc = oo the two terms enter the cross section incoherently. When the brem gluon has the 
same helicity as leg 4 and is coUinear with p4, it and gluon 4 are distinguished only by their p+ values. Then 
we arbitrarily call the one with smaller |p+ 1 the brem gluon. 
Now it is easy to see that 

Kk,4 = -Ptk {k+pif = -ptk^/k+ = -2ptk^k''/k+ (21) 



Then we have immediately, 

h 



dp4 dk ^1 ^ 2 I I A/\\2\ 



2b+|2|fc+|(27r)3' 



dP dk pt fik++p+)\ 2 ...X 

I^Corel (^^) 



2\P+\ |fc+|(27r)3 k++p+ \ p+2 {k++p+)y k 

The collinear divergence is now transparent in the integration over k near zero. It is the coefficient of the 
phase space factor dP/2|P+| that we should compare to the square of the tree amplitudes with self-energy 
corrections on external lines. 

We now show that the collinear divergence just isolated, when summed over all possible fc+ is canceled 
by a corresponding divergence in the self-energy correction to leg 4. This cancellation is a consequence of 
the Lee-Nauenberg theorem [7] , which stipulates that all collinear states with the same energy be included. 
The total energy of the brem gluon and the gluon represented by leg 4 is in the collinear limit 

fc- +P4- = ^ + A = A(l + kVpt) = (23) 
2k+ 2pJ 2pJ 2{k++pJ) 

which makes it clear that we should include all k~^,p^ consistent with fixed P"*" = k~^ +p^. We temporarily 
regulate the divergence by giving the collinear gluons a small mass /j.. Since lightcone phase space is mass 
independent, fi only appears in the above analysis in the expression for {k +P4)'^: 

(. + «f^-|t(i%M^S^) (24) 
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With this regulator the transverse momentuin integral we need is just 

/ =7rln ,,' „ (25) 

/o<fe^b+|<|fe+|A2 [k +^2(/,+ +p+)2/p+2]2 \k++pl\^n^e 

Then the coefficient of the jet phase space factor is 

The blowup as /i ^ is the collincar divergence wc arc seeking to resolve. According to the Lec-Nauenbcrg 
theorem, to get an infrared safe quantity we must sum over all in the range < \k'^\ < \P^\- And we 
must also include brem emission from the left of leg 4. The first term represents the emission of a brem 
gluon with identical hclicity to leg 4, so when we sum that term over the whole range of fc"*" we have included 
emission from both the left and right of leg 4. However the second term, represents the emission of a brem 
gluon with opposite helicity, and when summed over the whole range gives only brem gluon emission from 
the right of leg 4. The emission of an opposite hclicity gluon (with momentum k) from the left has the same 
squared amplitude, but it is convenient to switch the roles of k and p4, so fc always refers to the right gluon. 
Then the total emission rate is given by 



0<|fc+|<|P+| 



( \P^\ , |P+-fc+l^ , |fc+P V , , |fc+(p+-fc+)|A^ 

fiiv'^ \\h+(P+ -h+W^ li-+P+3| ^ \( P+ - h.+\P+3\ I ^ IP+|2,,2^ ^^'> 



A2|fc+|(27r)3' 

^ S>TT^\\k+iP+ -k+)\ ' |A;+P+3| ' |(P+-A;+)P+3|; \P+\''^?e 

Calling X = |fc+|/|P+|, |P+| times the quantity in parentheses can be rearranged 

1 (l-a;)3 a;3 / , x 1 

+ + = 2 ,t(1 -x) + 



x{l — x) X 1 — X \ 1 — X X 



f \k+\\p+ -k+\ lfc+l lfc+l 

~ \ p+2 + |p+_^+| + |p+_^+| ) (28) 



So with this notation the squared amplitude for jet production along gluon 4 is 



E 

0<|fe+|<|P+| 



4t^2 |p+| V ^ ' 1-x X J iJ?e ^ ' 

We still have to include the self energy corrections on the external lines. 
In Ref. [4] we obtained 

for the gluon self-energy for an off-shell gluon of momentum p after subtraction of the counterterms necessary 
to keep the gluon massless. Redoing the calculation with a mass /x for the particles circulating in the loop 
yields the modification 

^'^^ = |^^&' + -(i-->')(i+(T^ + i)i-((^' + -(i--)^''^^^^^ (31) 

- Constant +|^|l^(x(l- a:) + ^ + i^)ln(M^5e-+i) (32) 
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in the on-slicll limit 0. The constant (which is of order /i^) must be absorbed in a mass counterterm 

to keep the gluon massless. On the external line the is canceled by the extra gluon propagator, and the 
effect of the correction is just the wave function renormalization 

Of course the tree amplitude is multiplied by ^/Z for each leg and the squared tree amplitude by a factor of 
Z. Thus the correction on leg 4 is just 

{Z - l)|Aco.e|^ = l^^pf^ E (-(1 - + + ^) ln(M^<5e-+^) (34) 

Combining this self energy correction with (29) gives for the complete jet production, identifying P = p, 

We see that the coUinear divergence problem, which in lightcone gauge is confined to the self energy insertions 
on external lines, is resolved provided we interpret scattering amplitudes in terms of jets. 

However, this is not the end of the story because wc still see UV divergences (In S) and IR divergences 
due to X near or 1. As explained in [4] the latter divergences are regulated by discretization of The UV 
divergence from this calculation is to be combined with the UV divergences from triangle, box, and internal 
line self-energy diagrams to give the appropriate scale dependent coupling. Part of the IR divergences here 
will be canceled by soft gluon bremsstrahlung which we discuss next. In particular, we must find that the 
dependence on the resolution A is finite. But there will be residual IR divergences that are to be canceled 
by IR divergences from the other one loop diagrams (box, triangle, etc.). 

Our discussion of jet production only included the diagram with brem gluon attached the external leg 
identified with the jet (for definiteness we chose leg 4), neglecting its interference with diagrams with the 
brem gluon attached to other lines. This approximation is only valid, however, when the resolution is smaller 
than the momentum k^P4/p^ of the gluon in the jet. We can stipulate that all of the momenta pi are of 
order 0(1), but even so for small enough A;+ it is essential that we include the other diagrams. In this case 
all components of the brem gluon momentum are small and we must combine gluon emission from all legs 
coherently. In the large Nc limit life is simpler because we only need to include coherent emission from 
neighboring lines as already discussed. 

For definiteness, focus first on the coherent emission of a gluon between legs 3 and 4, both of which we 
assume to have outgoing helicity. Then the emission amplitudes are 



-2gAcore 



k++pt ^k++pt Klk 



= -2gA. 



Core 



k+p+ {k+pif k+pt {k+p3f 



(36) 



(37) 



k+{k+ +pt) {k + p4)^ k+{k++p+) (fc + P3)^ 

In these formulas we have assumed that Acore is the same in both terms, which is approximately true when 
all components of k are small. The case where k coUinear with one of the external momenta but not small 
will not introduce errors, because in that case the interference between different terms is negligible. For this 
we need to insist that no two external legs arc coUinear. which wc do. The squared amplitudes are: 

,f {k++p+fKl, {k++ptfKl^ {k++pt){k++p+)Kk,3-K,, 
l2A;+2p+2(fc + p4)4 2fc+2p+2(/c+p3)^ k+^p+ p+ {k + ps)^ {k + Pif 



iVi2 



= VIA 
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\2fc+2(fc+ 



p^Kl, , pt'Kls 



2fc+2(fc+ + p+)2(fc + P4)4 2fc+2(fc+ + pt)^{k+P3)^ 



k+\k+ + pt){k+ + pt){k+psf{k+pif . 
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These expressions contribute to the jet cross section for jets both along and p4. Let us define Vi = p^/pf . 
Then they should be integrated over the union of the two domains 1^4 : {k — k~^V4)'^ < k^A'^/p'^ and 
: {k — fc+vs)^ < k'^A'^/p^. In order not to double count we integrate over the entire first domain, but 
we integrate only over the part of the second domain satisfying (fc — ^+^4)^ > fc+A^/pJ. The necessary 
integrals over T>4 can be found in Appendix A. Define k = k — k~^V4, and give all of the gluons 3,4 and k a 
small mass /x. Then we have 

-2 
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2 {k+ + pjYii'^e 



(38) 



ptik-^+pJWe 



k'^Ptvl4 - A 



\k+\< 
\k+\> 



\pt\vl 



(39) 



k- {k — k~^vs4) 



In 
0, 



(fe + ii^k+ +pt)ypf){{k - k+v34r + i^Hk+ +pt)ypf) 

A4 



P4 ^34 



|fc+|> 



where the final forms are valid as /x 
{k + Pif = 



0. We have used the identity 

+ fi^ik++p+f _ p+\k-k+Vif+fl\k++p+f 



(40) 



(41) 



k+pf k+pf 

to simplify the integrands. 

In assembling these contributions we write separate equations for small and large k^, simplifying the 

coefficients in the first case: 

^'^ 9'\Acore\\j+'pfpfvi4{l-k+ptvl4/A^) „ A^ 



L 
L 



1V|2 



dk 



dk 



16|fc^ 

1V|2 , 



47r2|fc+| 



1^ 



A|2 



5' 



Core I 



16|A;+|7r3 



87r2 
/(fc++P3+)2 



(fc++pJ)2(fc++p^)VV 

{k++ptf , 



for \k+\ < 



\pt\v^ 



(42) 



34 



Pi 



pT 

pr 



pf 



{k++p+r 



(fc+- 

In ( 1 



-ptr 

A2 



In- 



fc+p+A2 



(fc++p+)Ve 



k+ptvl4 



for |fc+| > 



b4 l^^34 



87r2 



pf 



+ 



pf 



-21n 1 



{k++p+y 
A2 



In- 



fc+p|A2 

{k+TpfyJ^ 



for > 



A2 



\pi\v 



(43) 



34 



where the approximation in the last line is valid because the logarithm factor cuts off large fc"*". We still need 
to add the contribution of the part of the domain Vs that doesn't intersect 'D4. 

To handle the double constraint on the domain of integration, it is convenient to divide the contribution 
into two contributions, 



I : ^A^ <{k-k+V4)^ <'^Al; {k - k+vs)^ < '^A^ 
Pi Pi Pa 

II: ^Al<{k-k+V4f; {k-k+v^f <^A^ 

Pi pI 



(44) 

(45) 
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where we choose Aq ^ A large enough so that in region II we only need to include the diagram with the 
brem gluon attached to leg 3. It is not hard to show that in region II necessarily satisfies 




(46) 



Unfortunately the converse is not quite true: the condition on A;+ that implies k is in region II is slightly 
more strict: 



1 



''34 




(47) 



But at least the contribution from the part of region II that satisfies this last constraint is simply given; 



27rff^|Aco 



There remains the narrow window in k~^ 



In- 



k+p+A^ 



1 



^34 




{k+ +pJ)Ve 



< A;+ < ^ 

^^34 



for |A:+| > 



1 



+ 



A 



(48) 



^^34 



Ao 



+ 



(49) 



which contains a mixture of region I and II. But by taking Aq ^ A the contribution of this window can be 
made arbitrarily small, and at the same time A can be neglected in the lower limit on (48) 



87r2|fc+| V Pt^ 



+ 



pr 



{k++pt? 



In 



k+p+A"^ 



(fc++P3+)Ve 

for \k+\ > 



^0 



1^4^34 



(50) 



We limit the size of Ao so that all fc+ contributing to region / are negligible compared to the external p^'. 
In that case the integrand for region / simplifies to 



lV|2 



+ \A 



A|2 



VIA 



Core I 



-"■fc.4 



K 



fe.3 



Kk,3 ■ K 



kA 



k+^ik+pi)^ k+2(k+ps)^ k+'^{k+p^Y{k+p4Y 



(51) 



In region / one can show that -y/lfc+j > A(y|p^ — \J\pi\) / \j\p^l^\v:>A so we stipulate that > so 
that stays away from and it is safe to take the continuum limit of the k'^ sums. (To deal with the case 
\pt\ < 1^4 I ju^t switch the roles of legs 3 and 4 in the calculation). Now we use the identity (41) and 



to simplify the integrand even further 



k+^Pa ■ Pi - ^fc+p+(fc +P4)' - \k+pi{k+pif - fi^iptpt + k+pt + k+pt) (52) 



{fc+(fc 



Pt 



+ Pi) 



+ 



-ptk+jk+ps)^ - li^p+^ 
k+^{k+p3)^ 



2k+^P3 ■ Pi + k+p-^{k + pif + k+pt{k + p-if 
fc+2(/e+p3)2(/e+p4)2 

42,, ,2 f -M^^ • Pi \ 



(53) 
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Wc have dropped /i^ terms in the numerators when the denominators are prevented from vanishing strongly 
enough by virtue of being in region I. Since only {k + ps)^ is allowed to get small in region /, we only 
needed to keep when that factor appears squared in the denominator. The first term in braces will only 
receive contributions in the integration for fc^ = 0(/i^) as ^ 0. Thus the second constraint defining region 
/ collapses to a constraint on only: < k^p^v^^ < Aq, so the k integration can be freely done: 



27r / ' kdk-:^ = -tt + tt /' -tt = Trln -, for A^ < k+p+vl^ < A^ (54) 



The second term in braces is not only Lorentz invariant but only involves the variables constrained in 
defining region /; 

-2k-p3<A'^, A"^ < -2k ■ Pi < Al (55) 

(In this form the constraints coincide with the earlier ones only when fj, = 0. But the effect of the change is 
0(^1^) and negligible in region /.) Furthermore the measure for k integration dk~^dk/2\k~^\ = d'^kd{k^ + fx^) 
is also invariant. Thus it can be evaluated in any convenient frame (for instance, one in which P3 = P4 = 0.) 
The result is, assuming <C A, Aq <Cpi, 



•■j 

Jvi 



"^^k+psnk+p^y^ ^^"a^^"^^ ^^^> 



Of course we can't directly compare this to our previous results because they have not yet been integrated 
over A;+. But it is amusing to compare it with the integral of (50) over the missing range A^ < k'^p'^v^^ < Aq: 



f d\k+\ fc+ 
J \k+\ pt 



47r2 ./ |A;+| p+fj,^e A-k"^ A^ 



1 A' , . 
m — -r — - — h m ■ 



g'lAcoreW Al AoA3 



because v'^^ = (v^ — v^)"^ = —p'^i/p'^p^ ■ Remarkably, this integral exactly matches the entire effect from 
region /. The upshot is, that even though we didn't do the calculation this way, we can summarize the 
complete answer by quoting the following "results" for the fc integrals 

|^v|2 + |^A|2 ^ g2|Aco.cP [/ (fc++p+)^ ^ pr k+p+A^ 



16\k+\n3 8TT^\k+\ \ \ p+2 (fc++p+)2; (A;++p+)2^2g 

,(i^,^4L^),,^9E3^_„„(l_^)| ,„,n>^(59) 

V pt {k++P3rJ (fc++pj)2/i2e V k+pJviiJ } \PI\V34 

We must also remember that in executing the sum over k~^, the phase space measure is treated differently 
for the jet associated with each leg. Namely, the fc+ sum associated with the jet along leg i is taken 
holding A;+ + = fixed and there is an additional factor Ip^l/lfc"*" +pt\ = — A;+|/|Pj^| arising 
from transforming dp^dp^/p^ — {\P^ — k~^\/\P^\)dP^dPi/P^ . Fortunately these subtle modifications are 
only significant when A;+ = 0{pf), i.e. for a hard brem gluon whose contribution is dominated by a single 
diagram. Putting everything together we can write 

_ 5^ VlA^ 1^ V ' ( I iP^-k^f \,^ k^iPt-k+)A- 

^ 47r2 ^ \k+\ /x^e2 

\k+\<Ay\p+\vl^ 
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where we have used the cancellation 



47r2 



(61) 



valid after the (safe for these terms) limit of continuous A;+ . 

It is illuminating to rewrite (60) in a way that makes the symmetry under 3 4 manifest. 



i=3,4 |fc+| 



I gMcor.? sr^ 1 k+^vU\ptpt\ .... 

|/C+|<AV|P4+I"l4 



g^\Acorel_ 

^ 47r2 



2 



(63) 



Here we have picked m, the discretization unit, and have approximated 



E 



1 ^^ fc+2^t4|J-3+J'4+l ^ f^^/l^^"l'^-rft^^i2^4^|P3+p+| 



|/c+| A4 t A4 

A<|fe+|<AV|P+|,;32^ 

A2 A2 

= -In In — ^ (64) 

In the form (63) the symmetry 3 4 is transparent, but the finiteness of the A dependence, manifest in 
(60), has been obscured: the coefficient of In A2 has a small fc+ divergence on the first line that is canceled 
by a small k^ divergence on the second line. Another advantage of (63) is that the first line just gives the 
34 contribution to the production cross sections of jet 3 and jet 4 that matches our earlier discussion. In 
particular the fi ^ divergence is now transparently canceled by the wave function renormalization. 

A virtually identical calculation applies to the absorption of extra gluons in the initial state by the right 
of leg 1 and by the left of leg 2. Outgoing brcm gluons emitted between legs 1 and 2 are suppressed at 
Nc = 00 just as were incoming unseen gluons absorbed between legs 3 and 4 were suppressed. The result for 
soft gluon absorption between legs 1 and 2 is obtained from (60) by substituting 1,2 for 3,4. In this case 
k^,p'^,P2 are all positive so the many absolute value signs can be dropped. We note that this treatment 
of the initial state uses the Lee-Nauenberg procedure as a model of incoming legs as incoming jets, so the 
four legs of the core process are treated in a parallel fashion. This is in contrast to the by now standard 
procedure of absorbing the initial sate collinear divergences in the initial state parton distribution functions. 
This standard procedure is indeed approriate in interpreting collider experiments, where the gluonic process 
describes the scattering of the hard constituents of incoming hadrons. The Lee-Nauenberg procedure we 
follow is more general and works even in theories, such as A/" = 4 supersymmetric Yang-Mills, in which 
hadron-like bound states of constituents don't form. 

The situation for unseen gluon absorption and bremsstrahlung radiation on the left and right, either 
between legs 1 and 4 or between legs 2 and 3, is more complicated even at large Nc, because extra gluons 
in the initial state, the final state and both must be taken into account. This is because all these processes 
are now allowed at A^c = 00. For definiteness let's focus on the region between 1 and 4: gluons emitted or 
absorbed by the left of leg 4 and the left of leg 1. The diagrams for emission of a single unobserved gluon 
are shown on the second line of Fig. 1. The squared amplitudes are 

Pt^ (P4+-fc+)2] {k-k+Vi) 



|^V|2^^A|2 ^ 4g'\Acore\' 



2(P+ - fc+)2 2P4+2 



[(fe - k+V4)^ + t^^PtViP^ - k+f 
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Pi 



+2 



.2(P+-fc+)2 
P+{P+-k^ 



(p+ _ fc+)2 



2P, 



+2 



(fe - k+Vi^ 
[(fc - fc+Vi)2 + IJ?P+'^/{P+ - A;+)2]2 



(p+-fc+) p+{p+-k+ 

{k-k+vi)-{k-k+Vi) 



[{k - k+vif + /x2p+V(Pi+ - fc+)2][(fe - fc+'y4)2 + m'P4+V(^/ - / 

here we have used P^~ = pf + k~^, expressed the (fc + p,)^ in terms of Kij which we have written out 
explicitly. The main differences with the 34 contribution arc that the hclicitics of the two legs are opposite 
and Pi > while P4 , fc"^ < 0. Similarly with an extra unobserved soft gluon in the initial state the relevant 
diagrams have a gluon absorbed on the left of leg 1 or leg 4, but the result of the calculation is the same as 
for emission with the understanding that fc+ > 0. 

As we have already seen in the 34 case, the interference term is negligible unless all components of k are 
small so we can simplify that term by neglecting fc+ compared to the pf. For the future discussion we also 
write out separately the outgoing and incoming gluon cases sending fc — > — fc in the outgoing case so that 
fc+ > in both cases: 



I^Out 



2 , MA |2 
+ I^Outl 



VlAcorel' 
p+2 



P.' 



(P4+ + k+f 



+ 



2iP+ + k+f 



2Pt 



{k-k+vif 



[{k - k+Vif + ljfiP+^/{Pt + fc+)2]2 

(fe - k+vif 



2(P+ + fc+)2 2P+2 

2(fc - k+vi) ■ (fc - k+Vi) 



[(fc - fc+'i;i)2 + ^2][(fe _ k+ViY + /i2 

P4+ - k+f 



p+2 



P4+' 



+ 



[(fc - fc+i;i)2 + ljfiP+^/{P+ + fc+)2]2 

(fc-fc+i;4)^ 



(65) 



+ 



2(P+ - fc+)2 2P+2 
(P+ - k+f 



2P+2 



[(fc - k+Vif + li^PtViPt - k+ff 
(fc - k+vif 



2(P+ - fc+)2 

2(fc- fc+7;i) • (fc - fc+i;4) 

[(fc - fc+i;i)2 + /i2][(fc - k+v^Y + y?] 



[(fc - k+Vif + ljfiP^'^/{P^ - fc+)2]2 



(66) 



Clearly simply adding in these two contributions can't be the whole story since only one of them is needed 
to cancel IR divergences from loops. 

There is also a difficulty with the contributions for a hard unobserved gluon. A coUincar divergence is 
present when an extra hard outgoing gluon is collinear with either leg 4 or leg 1. When collinear with leg 4 
it is simply part of the jet associated with that leg and combines with the collinear gluon emission from the 
right of leg 4 to cancel the collinear divergence in the self energy correction to leg 4. The divergence coming 
from an outgoing extra gluon collinear with leg 1 has a very different meaning. Since we have stipulated 
that none of the gluons in the core process are collinear, this hard extra gluon is well separated from gluons 
3 and 4 and therefore in principle detectable: In this case the final state is unambiguously a three gluon 
state. A similar situation applies when an extra hard gluon in the initial state is collinear with leg 4. These 
non-jet-like collinear divergences have nothing to do with self-energy corrections on external lines and must 
be canceled by something else. 

The mechanism [7] that takes care of the doubled soft bremsstrahlung and the non-jet-like collinear 
divergences is shown in Fig. 2. At first glance it seems that these diagrams wouldn't be relevant, because 
they arc either disconnected (the first diagram) or apparently higher order (the diagrams involving two 
extra gluons). However, when we square the sum of these diagrams, the cross terms between the first term 
and the remaining three contribute as connected structures which are exactly the right order 0{g^) to be 
comparable to the one-loop and single gluon bremsstrahlung diagrams. The reason the last two diagrams are 
multiplied by 1/2 is explained in Fig. 3. It is essentially the same reasoning as that for multiplying self-energy 
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Figure 2: Diagrams representing an unseen bremsstrahlung gluon in both the initial and final state. The 
interference term, between the disconnected diagram on the top left and the remaining three diagrams, in 

the square of the sum of the four diagrams is of order 0{g^) and therefore comparable to the bremsstrahlung 
probability for one unseen gluon in the final state and none in the initial state or vice versa. To account 
for the factors of 1/2 see Fig. 3. 



bubbles on external lines by 1/2. Otherwise the squared amplitude would count equivalent contributions 
twice compared to what is required by closure for the unitary evolution operator U{ti,t2) or, equivalently, 
by unitarity of the S-matrix. 

We now discuss the evaluation of the cross terms in the square of the sum of diagrams in Fig. 2. Since the 
gluon momentum can't simultaneously be coUinear with both leg 1 and leg 4, the cross term with the second 
diagram on line 1 will only be singular for a soft forward scattered gluon. With incoming gluon polarization 
A this contribution is given by 



ip+ - k+){p+ 



l,k 



knpi - 



-kf 
{k- 



+ c.c. 



= -45- 



. jpt -k+)ip+ + k+)Ki^k-Ki,k 



[(fe - k+Vif + M2][(fe - k+vxY + M^] 



2 ~ k'^Vi) ■ (fc — k^Vi) 
^ {k-k+ViY{k-k+vxY 



(67) 



where the first form on the second line neglects fc+ compared to Pi and p^, and the second form also uses 
the fact that the IR divergence of this contribution is insensitive to the temporary gluon mass fx. Although 
the p+ dependent factors for the other polarization of incoming gluon V are slightly different, this difference 
disappears for soft gluons. So adding the two polarizations just multiplies this result by a factor 2. Comparing 
this expression to (65), (66), we see that it will cancel one of the two interference terms, so only one will be 
counted. 

Next we turn to the diagrams on the second line of Fig. 2. Since they are similar to each other, we only 
need do one in detail, say the first. It will give a singular contribution not only when k is soft, but also 
when it is hard and collinear with leg 1. In the latter case there will be three well separated gluons in the 
final state and therefore no confusion with the two gluon final state we want to describe. Nonetheless the 
collinear singularity from this contribution will cancel the one from a collinear brem gluon emission from leg 
1 with no extra gluon in the initial state. 

These diagrams with the forward scattering process entirely on an external leg are formally singular 
because the propagator between the last emission vertex and the rest of the diagram is on shell. In this 
regard it is analogous to self energy corrections on external lines. The most reliable way to handle such 
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-k \ 




k / 



Figure 3: One of the processes on the second hne of Fig. 2 viewed as a unitarity cut of a larger diagram. 
Because the unobserved gluon injects zero momentum into the single line it scatters from, the unitarity 
cuts A and B are equivalent and only one of them should be used in constructing the unitarity sum of the 
squared amplitude. This explains the factors of 1/2 in Fig. 2. The situation is entirely analogous to the 
well-known procedure of weighting self-energy bubbles on external lines with a factor of 1/2, as indicated 
on the right. In both cases the unitarity cut C is distinct and unique. 



situations is to compute the forward scattering process with pi off-shell, giving the intermediate gluons a 
small mass /i, and then go on shell by extracting the residue of the pole and factorizing it which means 
keeping only half of the correction to the residue (see Fig. 3). This process correctly discards the double 
pole contribution which is properly interpreted as an energy shift.'' 

Applying the Feynman rules to this diagram (without the factor of 1/2) with pi off-shell and k~ = 

''A rougher procedure is to attempt to work on shell from the beginning and use the ie's in the denominator to keep things 
finite. Then we would obtain 



,+2 



—ie[{pi + fc)2 + /i^ _ je] 
—ie[{pi + fc)2 + ij? — it] 



ICore 



(68) 



(69) 



When the cross term is constructed each of the above expressions will be multiplied by A^^^^ and added to its complex conjugate, 
and finally multiplied by 1/2: 



Cross Term^ 



Cross Term^ 



iVi^ii(p+ + fc+)2 



2 - 



1^2 fc+2p+2 I^Corel 



1 



1 



(pi + fc)2 -it (pi + fe)2 + /i2 + i 



Vlc2^AcoreP 



fc+2p+2 2[(pi 

^+2 



■fc)2+;^2_j^][(p^+fc)2+^2 



A;+2 



fc+2(p++A;+)2 p+2(p++fc+)2 



2[(pi + fc)2 + (u2_ie][(pi+fc)2 



(70) 



(71) 



Compared to the treatment in the text, which defines jet amplitudes as they would be extracted from larger diagrams with 
unitarity cuts, this rough procedure misses an overall factor (p]^ + fe+)/p^. For soft k this discrepancy is negligible, but for 
hard coUinear k the treatment in the text is the one that reflects the Lee-Nauenberg theorem. 
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(fe^ + iJ?)/2k+ we find 



(p+ + k+f 



Pi 



[plY[{pi + kY + ^j?] k+^pt^ 



Acorc{Pl,P2,P3,Pi) 



Ag\p+ 



-k+pt{k - k+vi)^ 



k+^pr{pl? 2[{k-k+Vi) 



2g\p+ + k+f 
k+p+ipir 



{k-k+vif 



(fe - k+vif + jj? 



- k 



+ 



^Core(Pl,P2,P3,P4) 



-{k++pt)pl/pt\ 

pI k+{k+ +p+){k-k+vif 



pT 



[(fe - fc+i;i)2 + /x2] 



212 



^Core 



(72) 



where the superscripts indicate the polarization of the unseen gluon. The first term in square brackets is 
a double pole in p\ with a coefficient that will not have a coUinear divergence (when smeared over fc.) Its 
residue will be proportional to the derivative of Acore but will not contribute to the collinear divergence. 
The residue of only the second term in square brackets, which is a single pole in pf will be divergent, and 
we easily read off the divergent contribution for this polarization 



2ff2(p++A;+)2fc+(fc++p+) {k-k+v^f 



k+p+ 



Pf 



[(fe-fc+t;i)2+^2]2 



^Core 



The other polarization is given by 
A'^ = -2g 



^k+{k+ 



Pt) 



pT 



pV 



k+^ 



k+{pt + k+f p^{pt + k^ 



{k~k+vif 



[{k - k+vif + fi^] 



r ^Core 



(73) 



(74) 



There are two contributions for VA polarization because the gluon connecting the absorption and emission 
vertices can have either polarization. 

When the cross term is constructed each of the above expressions will be multiplied by ^4^^^^,^. and added 
to its complex conjugate, which doubles it, and finally multiplied by 1/2 which undoubles it: 



Cross Term 



Cross Term 



AV 



^g'ipt 



(fe - k+Vif 



Pi 



[(fe. 



-2g' 



Pt 



+ 



k+Vif 



_pt + k+ ' pi^{pt + k+) 
Adding together the contribution for the two polarizations gives 

^{pt + k+f pt 



(fe - k+vif 
[(fe - k+vif + m2]2 



I^Co 



Cross Term 



-2g' 



Pf 



pt + k+ pnpt + k+) 



{k-k+vi) 



[(fe - k+vif + 



(75) 

el' (76) 

rl^Corel' (77) 



The process described here is one where an outgoing extra gluon is collinear with gluon 1 which is incoming. 
This gluon is thus not in jet 3 or jet 4. It can therefore be experimentally detected unless it is too soft. 
Comparing to (65) we see that the first two terms of (77) almost cancel the second term of (65). Integrating 
the difference over a neighborhood of fe = k+vi in the limit /U — > involves 



dk 



{k-k+vif 



[k-k+vif 



[(fe - k+Vif + ^2p+2/(p+ + fc+)2]2 [(fe _ k+Vif 



,212 



TT In ■ 



p+2 



(P+ + k+f 



(78) 



This vanishes as A;"*" ^ 0, when the extra gluon is undetectable, and when k+ is finite it will be excluded from 
an outgoing jet along or p4. Thus it shouldn't be included in the undetected bremsstrahlung associated 
with the core process. The last diagram of Fig. 2 removes the divergent contribution of the first term of (66) 
in a similar way. The last term of (77) is negligible for soft bremsstrahlung and, when k+ = 0(1) cancels 
a spurious collinear divergence from a brem gluon emitted from gluon 1 with the opposite helicity to that 
considered here. The squared amplitude for that opposite helicity process is 



I^Breml 



25^ 



(fe - k+Vlf\Acore{p\ - k,P2,P3,P4W 

P^iPi - k+f [(fe - A;+t;i)2 + ^2(p+ _ fc+)2/p+2]2 



(79) 
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Here the superscript denotes the polarization of the brern gluon, and —k is the incoming momentum of the 
brem gluon (so > 0). Notice that the first argument of ^core is Pi = Pi — k, not pi as in (77). To see all 
these cancellations, it is important to recall that k is to be smeared in a narrow region about the coUinear 
point k = k'^Pi/p^. Then P\ = p^ — k'^ {pf — k~^)pi/pf so the integration measure of the bremsstrahlung 
probability is 



dk dpi dk dPi P+ + k+ 



2\k+\2p+ ^ 2\k+\2P+ P+ 



(80) 



In (77) the first argument of ^coro is pi which is to be identified with p here. So we should compare (77) to 
(p'^ + k'^)/p'^ times the appropriate term in (65). 

After these cancellations what remains of the bremsstrahlung contributions (65), (66) is just 



2(fe 



:>+2 



+ 



(P/+fc+)^ 



2(P+ + fc+)2 
k~^Vi) ■ (k — k'^V4) 



[{k - k+vif + ij?][{k - k+Vi)"^ + At2] 



p+2 



2(p+ - k+y 



+ 



(p+ - k+y 

2P+2 



(fc - k+Vjf 

[{k - k+Vif + li^Pt^/{Pt + k+ff 



{k-k+vif 



[{k - k+viY + ^i^P+^/{P+ - k+YY 



(81) 



But now notice that these residual terms are the same as the 34 contribution if we identify P3 with —Pi. 
Thus when summed over k+ the 14 total contribution can be obtained from the 34 contribution with the 
substitution p3 — > —pi- This is consistent because unlike pi + pa, which is space-like, j>4 — pi is time-like. 

To summarize this section we collect together all the contributions from hard coUinear gluons, soft gluons 
and self energy corrections on external lines. The soft contributions boil down to a contribution like (60) for 
each pair of neighboring lines: 
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{Vi - Vjf = -{pi+pjf /P+P+ = \{Vi+P3f\l\PtPj'\- Notice that the terms on the first 
line correspond to contributions associated with each leg of the diagram, whereas those on the second line 
involve contributions associated with pairs of consecutive lines. The first category of terms includes the wave 
function renormalization due to self energy bubbles on external lines 
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(82) 



so that the coUinear divergence as ^ cancels. So we see that the temporary cutoff /x can be removed as 
soon as we combine everything together. We shall see that the remaining divergences in these expressions 
cancel against similar ones that come from the remaining one loop corrections to the glue-glue scattering 
process. These include self energy insertions on internal lines together with triangle and box diagrams. 

6 Cubic Vertex Function 



We shall not include calculational details for the one loop corrections to the cubic vertex function. They can 
be found in [12]. Instead we present the final answers for the vertex corrections with two on-shell gluons. 
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Figure 4: The triangle diagrams contributing to r^^^. The labels q,ko,ki,k2 are dual momenta. The 
actual momentum of any line is the difference of the momenta of the regions it bounds, pi = ki — ko, etc. 




We put the combination of swordfish and triangle diagrams (see Figs. 4,5) with two like-helicities and two 
legs on-shell in the form 



Tl loop + Tc.T. — 



■9' r 111 2 



(83) 



where the vectors ki,K carry the polarization of the two like- hclicity gluons, po is the four- momentum of 
the off-shell gluon, a = 1 when the on-shell gluons have like-helicity, and a = otherwise. Finally S is 
an infrared sensitive term that depends on the location of the off-shell gluon, but not on any of the gluon 
hclicitics. In the case ,P2 > 0, we denote by Sf (pi,P2) the value of S when leg /' is off-shell, and with 
loop momentum chosen so that is the longitudinal momentum of the internal line joining leg 1 to leg 3, 
satisfying < g+ < p^2- Then, 
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In addition to these corrections to the tree level cubic vertex, the triangle diagram with three like-helicities 
(see Fig. 6) is non-zero, and it is given, for the case of two on-shell legs, by 
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where Po is the momentum of the off-shell gluon. 



7 Box Diagrams 

For box diagrams the presence of q^ pole and double pole singularities in the integrand makes a direct 
evaluation of the integrals horrendous. Fortunately, it is possible to manipulate these integrands so that all 
of these problematic singularities reside in triangle-like diagrams. This is because on-shell tree amplitudes 
do not possess these singularities. We can identify tree amplitudes as sub-diagrams of one loop diagrams, 
but some of the legs of these sub-diagrams will be off-shell, so it would seem that features of the on-shell 



19 




limit can't be exploited. However, if one leaves the denominators of the off-shell tree subdiagrams alone, 

then the numerators can always be written as the on-shcll expression (with no q'^ singularities) plus terms 
each of which contain at least one factor of the virtuality g? of one of the off-shell legs. In a box diagram 
such terms will cancel a propagator reducing the required loop integrand to one with the structure of a 
triangle diagram. Since triangle integrals with (j+ singularities arc considerably easier to analyze than such 
box integrals, the resulting simplification is very useful. In the following subsection, we apply this technique 
to all of the helicity conserving box diagrams. (The helicity violating case was done in [4], where each box 
integrand was completely reduced to a sum of triangle-like integrands.) In the remaining subsections we 
complete the evaluation of the box diagrams. 

7.1 Box Reduction 

The thirteen box diagrams, seven for the helicity patterns A V AV, and six for A A VV, are shown in Figure 
7 and Figure 8 respectively. The integrand of any of these box diagrams has the structure 
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(89) 



where A/" is a quartic monomial of Kij^s carrying the gluon polarization information, and iZ is a rational 

function of q'^ ,pf . There arc only six possible A/''s: 

1. K^^K^5K^^K(>^: First diagram of Fig. 7. 

2. K^-i^K^r,K^^K^^: Second diagram of Fig. 7; first two diagrams of Fig. 8. 

3. K^-i^K^^K^^K^^: Third diagram of Fig. 7; third diagram of Fig. 8. 

4. K^-i^K^^K^^K^^: Sixth diagram of Fig. 7; sixth diagram of Fig. 8. 

5. K^^K^^K^^K^^: Fourth diagram of Fig. 7; fourth diagram of Fig. 8. 

6. K^^K^^K^^K^^: Fifth and seventh diagrams of Fig. 7; fifth diagram of Fig. 8. 

Since these structures come in complex conjugate pairs there are really only three essentially different struc- 
tures. The forms of the rational functions R change from diagram to diagram and we just list them in 
order: 

+2 +2 +2 +2 +2 +2 +2 +2 +2 +2 
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Figure 7: The boxes for the A V AV scattering process. The dashed boxes enclose sub diagrams whose 
replacement would convert the box to triangle-like loop integrals. 
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We can expand each of these thirteen rational functions in partial fractions and each will then be expressed 
as a sum of pure double poles, pure single poles, and a constant. 
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where is one of the four values 0,p'^ ,p^27 ^pt ■ Except for the fourth diagram of Fig. 7, one or more of 
the pole terms will be absent. Also C = 1 for the fifth and seventh diagrams of Fig. 7 and for the first and 
second diagrams of Fig. 8; and C = otherwise. 

The eight box diagrams with a hclicity violating subdiagram (enclosed by dashed boxes in the figures) 
can be completely reduced to triangle-like diagrams without coUinear divergences. Each of those completely 
reducible diagrams has one of the first four polarization structures in our list, i.e. two neighboring K's have 
the same polarization. Then for a like-polarization pair, one can use an identity like 
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Figure 8: The boxes for the A A VV scattering process. The dashed boxes enclose sub-diagrams whose 
replacement would convert the box to triangle-like loop integrals. 



which underhes the on-shell vanishing of the three like-helicity amplitude to convert the integrand to a 
triangle-Hke one which is free of collinear divergences. The integrals over q,q~ can be evaluated as in [4]. 
The sum over can be converted to an integral and carried out for the terms with no poles in . The 
sums with g+ poles are left undone as what we call the infrared divergent contribution, and will eventually 
be canceled against real gluon bremsstrahlung in their contribution to cross sections. For more detail on the 
calculation of these eight diagrams, sec Appendix E. 

The remaining five diagrams cannot be completely reduced to triangle-like diagrams. But we can ma- 
nipulate the integrand so that all g'+ divergences are located in triangle-like diagrams. Then the remaining 
box integral can be straightforwardly evaluated. Since the procedure works in essentially the same way for 
each of these five diagrams, we shall illustrate the method by picking one of the two polarization structures, 
the last in our list, which appears in the fifth and seventh box diagram of Fig. 7 and in the fifth diagram 
of Fig. 8. As already mentioned the constant term in R only appears in these two (the fifth and seventh of 
Fig. 7) of the five "difficult" diagrams. Its contribution is the same for both diagrams, the integrand being 

1 K'^ R'^ 
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(27r)4 p+p+p+pt{q - konq - hfiq - k^fiq - h) 

and the integral of (93) is evaluated in Appendix D. 

Of the pole terms in R it is sufficient for our illustrative evaluation to single out only one of the q~^ pole 
locations, say that at g+ = pt ■ Thus we wish to integrate the following integrand 



(27r)4 ^ (q+ - p+y q+ -pt J ptptptptlhUhl 



(94) 



where we have introduced the shorthand notation gj = q—k^, with ki the dual momenta related to the actual 
momenta by pi = fcj+i — fcj. We usually take k^ = but leave feo arbitrary. Next we list eight identities, 
the last two of which enable the desired manipulation of this model integrand. 

K^,K^, + K^,K^, = K^, . 7^64 = + f (P^ (9+ + pt) + ptipt - 9+)) 
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K^,K^,,-K^,Kl = ^(iff^i^e^'i-^iWi+^iW-^iWJ (96) 
Pu 

I (g+-Pl2)P2"g3 _ (g+ - Pl2)P3 g? (97) 

^^5^2^5-^2^5^3''5 = ^i^^ (if3'^2^5''2 " ^3'2^5^2 + ^3'^2^5'3 " ^3^2^5'^3) (98) 

P23 

K-,K-,+K^,K^, = K,,-K,,= ^P^'^''p'P'' -'I{ptq++pt{pt2-Q^)) 

I (g+ +P4)P^gg _ (g+ +P4)p|gi (gg) 

J^35<4^<4^3"5 = ^^^{K,\K-,-K^,K^,) (100) 
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^ iPl - g+)P2 go _ (Pi^ - g+)pi^g2 (^Q^) 

= ^^^(ii-2W5-i^2W5) (102) 
P12 

The first six identities are needed for other box integrands. All but the first terms of any of the right 
sides contain a factor of qf which would cancel one of the propagators converting the box to a triangle-like 
diagram. Depending on which term in the partial fraction expansion we consider, we can use one of these 
identities to switch A and V on a pair of K factors. For our model case we use the last two identities to 
write the first two K factors in two diff'erent ways 

^aW5 = -K^.K^. + - f + |(P2^g^ +Ptipt. - g^)) + - f - -f 
K^,K^, = K^,K^,+ ^—^{K^,K^,-K^,K^,) (103) 

Pl2 

We use the first rewrite for the double pole and the second for the single pole. The first terms of either 
convert the polarization structure to one which completely reduces to triangle-like diagrams free of collinear 
divcrgcncx's just as with the eight box diagrams with hclicity violating subdiagrams. Their calculation follows 
the models given in Appendix E. The second terms of either have an explicit factor that cancels the singular 
g+ denominators, leaving a box integrand free of g+ divergences. However, only the box from the second line 
is free of collinear divergences; its evaluation is therefore straightforward and is given in Appendix F. There 
arc no more terms for the second rewrite, and the remaining terms of the first each have a virtuality factor 
go'gi'gl which converts the integrand to a triangle-like one. Unfortunately, both these last triangle-like 
diagrams and the box integrand left after the first rewrite individually have collinear divergences. They 
must, of course, cancel among themselves, but to achieve clean results we discuss in the next subsection a 
way to rearrange the integrands to finesse this diSiculty. 

Summarizing this subsection, we have rewritten our model integrand in the form 
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In the next subsection we show how to deal with the coUinear divergence in the first term of the right side 
and in the triangle-hkc diagrams associated with it: 
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(105) 



ActuaUy the prefactor kills the coUinear divergence in the first term in the square brackets so we only need 
deal with the last two terms. 



7.2 Subtracting CoUinear Divergences 

The box reduction we have so far accomplished has the undesirable feature that the new box integrands 
have coUinear divergences that were not present in the original box integrands. This means that there 
must be canceling coUinear divergences among the triangle-like diagrams that we generated in the reduction 
procedure. The terms with four K^s in the numerator will not have this problem but all the terms with only 
2 K^s in the numerator do. In order to deal with this problem we must add some triangle-like diagrams to 
these problematic box integrands in such a way as to regulate these divergences. To begin, we note that 
the terms linear in loop momentum dependent K^s can be made IR finite by a simple subtraction of two 
triangle-like terms. We find that each of the combinations 
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is finite integrated in the infra-red, and the triangle-like subtractions are quadratically convergent in the 
ultra-violet. This nice IR behavior is not spoiled by multiplying each expression by further factors of loop 
momentum dependent K's, and up to two such factors could be applied keeping the UV behavior no worse 
than logarithmic. Thus we can satisfactorily regulate the terms quadratic in loop momenta by simply 
multiplying one of these expressions by the appropriate K. There are several choices for regulating each 
term, so we arbitrarily choose one of them. Of course, as already mentioned the term quartic in the loop 
momenta is IR convergent by itself and needs no subtractions. 

When we pass to the Schwinger parameterization (with the notation of Appendix C) the numerator 
factors in (106) enjoy a nice simplification after an appropriate shift in q: 
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Here Ko = X2P1 + .13 (pi +P2) ^ x^p^ is the 6^0 limit oi K — ko where K has been defined in Appendix A. 

Now consider the terms quadratic in loop momentum dependent K^s. These terms involve one of the pairs 
{K25,K35), {K35,Kq4), {Kq4,Kqi), {Kqi,K25), each pair is associated with a neighboring pair of external 
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lines of the box diagram. The two members of each pair have opposite polarization, with both possibilities 
occurring. It is easy to confirm that the terms involving each pair can be obtained from one another by 
cyclic symmetry. Therefore we need only analyze one class of terms, say, those involving the pair occurring 
in the model integrand of the previous subsection, K^^,Kq^. These terms can be regulated in the infra-red 
either by multiplying the last line of (106) by K^^ or the third line by Kq^. We choose the former and for 
the other pairs make the choice dictated by cyclic symmetry. The terms in the numerator that will survive 
integration over q are 

K^MpI,pI, - KX.pUI - K^,pWo) - -^q^'^ixrK^ + X2K^:,){K^^p\^ + K^^pl^) 

+{x2K^i + xiK^J{KXyi2 + K^pl4){xiX3pl2 + X2Xipli - q^) (108) 

In addition to dropping terms that directly integrate to 0, we also used —q~^q~ q^q^ valid under q 
integration. As shown in Appendix C, the effect of in the numerator is its replacement by a factor 
2H = 2(xiXzpi2 + X2XApii) and the effect of g^g^ is its replacement by H/2. Thus we can replace 

K^Jk^,,-^-^) -2{x,xsp12+X2X,pI,){x2K^,+x,K^,)(^ + ^) (109) 

V Pl2 Pl4 J \ Pl4 Pl2 J 

= -4(xiX3pL + X2X4p\a){x2K!^^ + x^K!^^) ^l^^^f\ (110) 

P\PtP\2PiA 

The loop integral involving this pair of ii''s can be done (see Appendix C) 

/■ pi2K^M r = / d'^q Pl2K^5 [K^A - Kt4qllp\2 - KlAlpj^ 
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The box integrands involving other pairs of K^s can be regulated and evaluated in an exactly parallel fashion. 
We shall quote the results of combining all contributions in our results Section 9. 




7.3 Calculation of the Triangle Diagrams with CollineEir Divergences 

We now turn to the triangle-like diagrams which contain coUinear divergences, the last two terms of 
(105). These divergences must be canceled when we add back in the triangle- like diagrams we subtracted 
from the box to cancel its collinear divergences. 
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The last two terms in the square brackets on the right have a collinear divergence due to the vanishing of q^ 
and qf when q is collinear with pi, specifically at g = q'^Pi/pt ■ In this limit we have 
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We see that this collinear divergence cancels. Similarly the first and last terms have a collinear divergence 
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and they also cancel. The other possible collinear divergences in these expressions are killed by the polar- 
ization factors ifg4 and if 35. 

To calculate these triangle-like contributions, we remove from each triangle structure its collinear limit. 
We work out one case, the first one in (113), explicitly. We integrate dq~ first*. The q' integrals of all six 
types of bubble integrands and all four types of triangle-like integrands are listed in Appendix B. Note that, 
under the assumption that pf < —P4, for < < pf only the qo pole contributes, for pf < q~^ < —pt 
only the 93 pole contributes, and the integral gives otherwise: 
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for <q~^ <pf 
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We see indeed that no collinear divergences are encountered when these expressions are integrated over q. 

The other cases are similar. The fact that the 012 integrand has two collinear divergences that require 
cancellation is not a problem, because, as seen in appendix B, after the integral over q~ , the two divergences 
occur in disjoint regions of q'^ , so one can simply remove them additively: 
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Another approach to the collinear divergence problem is to calculate the individual diagrams with a mass 
regulator /i^, and send /i ^ only after combining the terms that together are free of collinear divergences. 
In fact, for the purposes of automating our calculations, we found this latter regulator method more efficient, 
and used it to obtain the combined final results of all these integrations quoted at length in Appendix G. 



8 Non-box One Loop Corrections to Scattering of Glue by Glue 
8.1 Cubic Vertex Corrections and Self Energy Insertions on Internal Lines. 

We quote here the contribution of triangle corrections to glue-glue scattering combined with the self energy 
bubbles on internal lines (see [4, 12]): 
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®A11 our formulas have treated d'^q as Euclidean, but to integrate over q we convert to Minkowski measure d'^q = —id'^qM, 



and use the usual prescription for propagator denominators 
sum of residues and the i's cancel. 



ie. Then the integral over q just gives 2m times the 
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9+ 

8.2 Qucirtic Triangle Diagrams 



-S2{-P4, -P23) - Si{pi4,P2) + >t(P?4.Pl4) 



(119) 
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Figure 9: The quartic triangle diagrams shown gonerically, without arrows indicating spin flow. Particle 
labels 1234 are applied counter-clockwise starting at the lower left of each diagram. 

There are four distinct quartic triangle structures (see Fig. 9), which we label by the two legs entering 
the quartic vertex. Half of the diagrams for each polarization configuration arc given in Figs. 10, 11. The 
integrand of each diagram has three of the four possible propagator factors 1/qf for z = 0, 1, 2, 3, a numerator 
consisting of one of the eight polarization structures Kq^Kq^, Kq^Kq^, K^^Kq^, K^^Kq^, K^qK25, iirggif^g, 
A'l^A'^r,, K^^K 2r-,, times a rational function of and the The q~ and q integrations are virtually 
identical from one diagram to the other. 

A model integrand for a quartic triangle is: 

-K- K- 111 

1 4 

^0+ < a+ < kt \[\ogH, + log(5e-)]^^^^|^^^p^ - i [log if, + log(5e-)](fc3+ - g+) 

kt < g+ < kt \ [log + log {Se'')] (-^3+ + 4)(fc| - g+) ^ 1 [^^g ^ (^^.^j _ 
4 {k^ - kl) 4 



kt<q+< kt ^ [log + log {de->')] ( ^3+ + fc2+)(fco+ g+) ^ 1 [i„g ^ i„g (^g.)] _ 

(Kg — j * 
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Figure 10: Half of the quartic triangle diagrams for the AVAV scattering process. The six others are similar 
but with the quartic vertex at the left or top. Particle labels 1234 are applied counter-clockwise starting 
at the lower left of each diagram. 



The various H's are defined: 



k+ <q+ < k+ 



kt <Q~^ < ki 



k^ < < k^ 
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Hd = 
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Hs = 
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pUkt-q^){kt-q^) 


{ko — k2){^k^ + 
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{—kt + kt)(kt — kt) 


pUk+-q+){-q+ 
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{kt — ko){kQ — 


kt) 


{kt ~ kt){kt ~ kt) 


P?2(9+-fco+)(4- 


9+) 


Hd = 


Pl2{kt-q+){kt-q+) 


{k+-k+r 




{—kt + kt){kt — kt) 


Piiikt - <7+ 


? 


Hu 


pU-q"- +kt){kt -q-^) 


(— + ki){—k'^ 


+ kt) 


{—kt + kt){kt — kt) 



Assume that the coefficient of this diagram is A + B{q'^ — fcjl") + single pole + double pole. Multiply the above 
result by the prefactor, partial fraction the coefficient of the logarithm into pole terms and polynomials. The 
polynomials can be integrated to give:: 

^{kt - kt){kt - kt){-4Bkt + 2Bkt + 9A + 2Bkt) 
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Figure 11: Half of the quaxtic triangle diagrams for the A A VV scattering process. The five others axe 
similar but with the quartic vertex at the left or top. Particle labels 1234 are applied counter-clockwise 
starting at the lower left of each diagram. 



^^{Bk+ - 2Bk+ + 3A + Bk+){kt - k+){k+ - k+) log ipj^Se'') 

It's quite interesting that although in principle, a single pole term in the prefactor can still contribute, their 
net effect is zero. After the polynomial terms are integrated out, the remaining terms will be the infrared 
terms, they either cancel or combine with infrared terms from other diagrams into complete trees. 
Similarly, the quartic triangle integrand 

111 



^1,6^5,2 



if Qo ll 



gives: 





-q+){kt- 














k^) 




{k2 — k^) 




{kt 







k+<q+<k+ + -[logH, + \og{6e'') 

kt <q^< kt + \[logH, + log (fc^)] ^"^^ '^^ ' - ^[logiJ„ + log ((5e^)](fc+ - q+) 

kt<q+<k+ - l[logg. + log (fe-^)] "(fci^i'fcl)^ + i[logg^ + log i^^'mt 

Assume that the prefactor of this integrand is A + B (5+ — ki)+ single pole + double pole, a similar procedure 
gives for the A, B terms: 

^(-4 + kt){-k+ + k+){-4Bk+ + 2Bk+ + 2Bk+ + 9A) 

+ ^{kt - kt){kt - kt){Bk+ - 2Bk+ + Bk+ + 3A) XogHpl^Se^)) 

Next, the integrand 

-K- 111 
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gives: 



k+<q+<k+ + -[\ogHr + \ogide')]{k+ -q+) 



<q^ < k^ 







Assume the prefactor of this model integrand is: A + B{q~^ — k^) + single pole + double pole. Then the A, 
B terms can be integrated to give: 
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{kt - kt){-k^ + kt)i2Bkt + 2Bk^ + 9A- ABk^) 



Finally, the integrand 



+ ^(fco+ - ktKkt - k^KBk+ - 2Bk+ +3A + Bk+) log {plje'') 



-K- K- 111 



gives 



k-^ <q+ <k^ 

kt<q-^<kt + \[\ogHt + log (fe^)] ^-^^ ^fcf!!~f ^ - \[^ogHi + log - g+) 

k+<q+<k+ - -^[\ogHi+\og{5e^)]{k+ -q+) 

Assume the prefactor of this model integrand is: A + B[q^ — ^2^) + single pole + double pole. And the A,B 
terms can be integrated to give: 



l(fc+ - k+){kt - kt){2Bk+ + 2Bk+ - Ak+B + 9 A) 
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{kt - kt){kt - kt){Bkt + Bkt + - 2fc+B) log {pl^Se^) 



8.3 Double Queirtic Diagrams 



A typical double quartic integrand is: ^ \ This gives: 

9i 93 



87r2- 



kt 



dq+ 1 

o log 



(9+-fc+)(A;+-g+) „ 



{kt-k+Y 



We can also have a diagram like Its contribution is similar, but spans over all three regions. 



dq+ 1 

kt ^2 — kg 2 



log 
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Figure 12: The double quartic diagrams for the two possible polarization patterns. Particle labels 1234 are 
applied counter-clockwise starting at the lower left of each diagram. 



Assuming that the prefactors of the double quartic diagrams are A + pole terms. We can integrating out 
the A term, leaving the rest as infrared terms. Thus we have: 

_i__l[_2 + log(p2^fc'^)] 

for the first case and 

_L_l[_2 + log(p?25e^)] 

for the second case. 

9 Final Results 

In the previous sections we have described our calculational methods by choosing a single example of each 
distinct type, and analyzing it in detail (relegating the more tedious parts to appendices). These examples 
are chosen to illustrate every type of technical complication we encountered. However, along with each 
such example there are quite a few others involving essentially identical calculations. In fact, there are so 
many that we chose to automate their calculation using Matlab and Mathematica. After the results of all 
these many calculations are combined, there ensues a stunning simplification that allows us to present the 
complete elastic glue-glue scattering amplitude in the first subsection, and in the second subsection, to obtain 
the complete answer for probabilities including unseen gluons in the initial and final states. This last result 
is compact, infrared finite, manifestly Lorentz invariant, and displays the ultraviolet behavior dictated by 
asymptotic freedom. 
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9.1 One Loop Corrections to Elastic Scattering of Glue by Glue 

We quote here the amputated four gluon amplitudes, which do not include any external leg corrections. 



A 1 — loop 
^AA VV 
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The infrared terms for both helicity arrangements are the same multiples of the corresponding trees: 
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The reader will note that the infrared sensitive terms depend on the ultraviolet cutoff S. This entangling of 
infrared and ultraviolet divergences is a familiar consequence of the way we have cut off p+ = singularities. 
These entangled divergences are precisely cancelled by similar divergences in the self-energy corrections to 
external lines which contribute to the Y[ V^i factors that convert the amputated amplitudes to properly 
normalized scattering amplitudes. When these factors are included (as they will be in the next subsection on 
physical probabilities), the net coefficient of In 5 becomes — llg-^/247r^ in precise agreement with asymptotic 
freedom [8,9]. 

The terms in these amplitudes that are not multiplied by trees are Lorentz violating anomalies that 
must be removed by counterterms. In Section 10 we show how this can be done locally in target space 
and described locally on the worldshcct. As wo shall see, after these counterterms are taken into account, 
the only change in the rest of the formula is a change 73/9 — » 67/9 in the constant terms multiplying the 
respective trees. We assume these changes have been done in our discussion of unseen gluons in the following 
subsection. 

The expressions for the loop amplitudes are real in the unphysical region for scattering in the 12 channel 
where = —s and = —t are both positive. The physical region for this process is s > 0, t < 0, which 
we can obtain by analytic continuation. Since the physical s is above the cut on the positive real axis, we 
obtain the physical amplitudes by substituting = s-^^^s in the above formulas. In this way we see that in 
the physical region the amplitudes acquire an imaginary part. In the next subsection we use these physical 
region amplitudes in the calculation of gluon detection probabilities. 



9.2 Probabilities Including Bremsstrahlung and Unseen Initial Gluons 

In this section, we will focus on the case when the extra gluon(we can break the Bose symmetry by 
defining it to be the softest one) is between particle 3 and 4. Recalling the results of section 5, there is a 
total of three terms that contribute to the infrared and coUinear singularity. 
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While the infrared terms from loop calculation in region 34 is: 
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First, notice that Mcoii is organized according to which leg the colhnear emission is attached, but for the 
sake of the current discussion, we need to break them up into different regions. For example, rewrite the 
term (take i = 3) 



in+l 



fc+(|P3+| 



k+) 



k+\p+\^ 



(\p+\-k+)\p+\^ 



as 



1 



k+f 



+ 



+ 



\k+Pt\ \P^\-k+ {\P^\-k+)\P+'\ 



The first two term are related to the last two by substituting fc"*" {P^ — The first two will be 

divergent when the extra gluon becomes soft. The last two will be divergent when the extra gluon becomes 
dominating over gluon 3, so, by the definition of 'the extra gluon' given above, gluon 3 becomes the 'extra 
one'. Hence we assign the last two terms to region 23. 

Second combine Mioop with Mgoft brem (leaving out a common factor ^ |^core| ) 
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We can see that the first two terms in Mgofi f,rem cancel the divergence in the first two terms of Mi^op- 
Performing the first two integrals: 
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Rewrite the divergent integral as 
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The first two integrals will be cancelled by Mcoii later. We get: 
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Putting everything together, we get: 
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Thus, after including the (12), (23), and (41) cases, the total probabilities are: 
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We see that all IR divergences have cancelled, that the UV divergences are exactly those dictated by asymp- 
totic freedom, and Lorentz invariance is manifest. 



10 Worldsheet Description of Counterterms 

As we have seen, counterterms that are polynomials in the dual momenta must be added to the two, three 
and four point functions in order achieve the correct Lorentz covariant results. Specifically in [4] we required 
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the following counterterms: 
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which are all polynomials in the dual momenta. 

Quartic counterterms must also be included in the list. At first sight the terms in F that need to be 
canceled seem to be rational functions of the pi (which would be nonlocal in a;~): 
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Local counterterms must be polynomials in the momenta. However, we note that the nonpolynomial parts 
of these anomalies are proportional to the quartic vertex contributions to the corresponding tree amplitudes. 
The addition of a term to n^Yp proportional to p^, which is an allowed counterterm by power counting, would 
contribute a term proportional to the part of the same trees built from pairs of cubic vertices. By tuning 
the coefficient of p^ we can convert these nonpolynomial anomalies to complete trees, obviating the need to 
cancel them with a counterterm. They just correspond to a perfectly allowed finite coupling renormalization. 
The change in II^Yp which accomplishes this is just 
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After this rearrangement we see that the required quartic counterterms are simply constants: 
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It is noteworthy that these quartic counterterms are spin independent. This is consistent with the inter- 
pretation of the anomalies as ultraviolet artifacts of box diagrams: The large momentum behavior of box 
integrands must of necessity be of the form qi^^ qi^^ qi^^ qt^'^^ J which integrates to completely symmetrized 
Kronecker deltas. 

This is all quite satisfactory from the field theoretic point of view, but the worldsheet description makes 
more stringent requirements on the counterterms: They must be generated by purely local changes to the 
worldsheet action. Worldsheet locality is quite independent of field theoretic (target space) locality, and we 
must still show how it can be preserved. We start with IIq^ . Interpreted as a contribution to the worldsheet 
path integral representation of a gluon propagator, it should be multiplied by T'/2p+ [4]. 



2p+ 



167r2 



/ 



dr 



g^2(0)+gA2(p+) ^2 



P-" 



487r2 



J drda ^ 



(130) 



When exponentiated (through higher loop corrections) this expression can be interpreted as adding new 
terms to —S, where S is the worldsheet action. The first term modifies the treatment of the worldsheet 
boundary in a similar way to the description of a mass term [5], and the second term is a new bulk term. 
Both terms violate helicity in the right way to cancel the helicity violation implied by the nonvanishing of 
n^^. These new terms also produce new effects from contact contributions arising when the bulk term sits 
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on the same time; slice; as an interaction vertex, similarly to the generation of quartic vertices from two cubics 
[3]. Using the generating function obtained in that reference we find, for correlators on the same time slice, 
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Here o, to are the discretization units of t, a respectively. The terms proportional to a will be negligible in 
the continuum limit, so we see that the only contact contribution that survives is the last term on the last 
line. Its effect is to produce new quartic vertices similar to the quartic anomalies already discussed. Indeed 
the usual quartic vertices arise from the second term of the correlator 
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which has an exactly analogous 1/a contribution. By a parallel calculation we easily find that the contraction 
terms produced by 11^^ give the following quartic vertices: 
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These new terms add to the anomalous contributions previously discussed: 
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The nonlocal terms in these expressions can be handled as before by retuning the p^ term in 11^^ a little 
differently: 
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The worldshcet description of the first two terms in n^Yp has been explained in [4] : the first term can be 
absorbed in a worldshcet bomidary cosmological constant, and the second term, which has the interpretation 
as a shift in the gluon (mass)^, can be absorbed in the worldshcet description of a mass counterterm [5]. 
The way to put the last term in the worldsheet description is a little more subtle. We first note that its 
effect is simply a finite contribution to wave function renormalization Z which multiplies the cubic vertex 
by and the quartic vertex by Z^. It is convenient to, at the same time, make a finite renormalization 
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of the coupling c;onstant g g/Z to reduce the effect to multiplying the cubic vertex by Z^^^ leaving 
the quartic vertex untouched. Then the net effect, to be described by the worldsheet formalism, is to 
modify the constant part of the correction multiplying the part of the tree involving gluon exchange from 
73/9 67/9 without touching the correction multiplying the quartic vertex part of the tree, which gets 
adjusted to 67/9 by the nonlocal part of the quartic anomaly. One's first thought is to simply change the 
coefficient of the cubic vertex appropriately. But then the worldsheet contact contributions would make a 
corresponding modification to the quartic vertex contribution and the change would only amount to a finite 
renormalization of g. Fortunately it is possible to prevent the dq/da worldsheet insertions from generating 
quartic contributions by altering the ghost worldsheet action near the interaction point. 

The crucial feature of the ghost path integral that permits this was explained in [2]. The discretized 
(p+ = Mm) ghost action on a fixed time slice is 

M-l 

Y,ibi+i-h){c^+i-Ci) (141) 

where bo,bM,co,CM = 0. The effect of the worldsheet integral of the exponential of this expression is to 
supply a factor of M which cancels a 1/M from the coordinate part of the path integral. If a single link in 
the sum is deleted, the result of integration is down by 1/M in other words it is 1. This was the mechanism 
we used to generate needed l/p"*" factors in the vertex functions. The location of these ghost deletions is 
indicated by short vertical lines in Fig. 13. If two links are deleted on the same time slice of the same gluon 



^1 



Figure 13: Discretized worldsheet for a four gluon tree. The solid squares indicate where dq/da insertions 
can be located. The short vertical lines indicate the links to be deleted in the worldsheet ghost action. All 
indicated ghost link deletions are present regardless of the insertion location. 

propagator, the worldsheet integral gives zero. Thus, when the insertions are on the same time slice there 
will be two deleted links and the contribution will be suppressed. Since the deleted links produce unwanted 
l/j3+ factors, one must also include dummy ghost insertions (defined in [3]) to put back corresponding 
factors of p+. Thus we can suppress contact contributions from being produced by the cubic counterterms 
by accompanying each dq/da insertion with an extra deleted link. We show in Fig. 14 which ghost link 
deletion is made for each of the six possible dq/da insertions (three for the fusion vertex and three for the 
fission vertex). This ghost deletion scheme allows us in eff'ect to change the cubic vertex by a multiple of 
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Figure 14: Discrctizcd worldshcct for a four gluon amplitude with one tree cubic vortex (solid squares) and 
one cubic counterterm vertex (open squares). The squares indicate where dq/da insertions can be located. 
The short vertical lines indicate the links to be deleted in the worldsheet ghost action. Notice that each 
counterterm insertion is accompanied by an extra ghost link deletion. Inspection shows that whenever two 
insertions are on the same time slice of the same gluon propagator, there are two deletions and hence the 
contact contribution is suppressed. 
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itself without affecting the quartic vertex, and this in turn allows the conversion of the nonlocal parts of the 
quartic counterterms to complete trees. 

The same scheme is very useful in translating the cubic counterterm (125) to the worldsheet. 

3 3 

^cCt. — + ki + = -|^2^2 [^0 ki + k2 — ki + 3fcf ] 

= :£5K-pf+3<l = ^(pj2^(B)-rf|^,^) + V(.4)) (142) 

Here A and B label worldsheet points just to the left and right of the internal boundary separating the two 
gluon propagators that fuse to or fission from the third gluon propagator (see Fig. 15). We can then use 
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Figure 15: Worldsheet for cubic fusion vertex. 



the ghost deletion scheme just described to guarantee that these insertions produce no modification of the 

quartic counterterms. 

We have now shown how the self-energy and the cubic counterterms, together with the nonlocal parts 
(in target space) of the quartic counterterms can all be given a local worldsheet description. It remains to 
find a local worldsheet description of the purely constant parts of the quartic counterterms 

pAAVV _ n -pAVAV _ ff^ f-IA'i\ 

^ C.T.ws — ^5 C.T.ws — y'-^'^) 

We cannot simply postulate a direct quartic interaction vertex in the worldsheet formalism without destroying 
worldsheet locality. We therefore search for a cubic vertex whose contact contributions generate constant 
quartic vertices. Consider the simple ansatz 

C-^v =5^e(P2" -Pf) - .9-^e(p^^(2) -P+^(1)) (144) 

where the legs of the vertex are labeled 1,2,3 counterclockwise and 1,2 have like helicity. Then it is not 
hard to see that the four gluon trees built from one tree cubic and one of these vertices, generate the contact 
contributions 

^AAVV ^ _gi^^ ^AVAV ^ (^45) 

That is the ratio of the two polarization structures is the same as that coming from the term, 
in the field theoretic Lagrangian. To make this work we need to suppress the new cubic couplings while 
retaining the contact contributions. To do this we can write 

= C^^^ - C^^^ (146) 
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and apply ghost link deletions on the second term so that it will not produce contact contributions. Then 
the exchange graphs will cancel leaving only the contact quartic vertex! 

We have found no simple variation of this scheme that provides us with exactly the counterterm (143). 
Instead, our proposal is to increase the flexibility of the worldshcct formalism by increasing the dimensionality 
of the worldsheet fields q{a, r). This is not unprecedented. Recall that in the AdS/CFT correspondence [16], 
the string theory is formulated in ten space-time dimensions whereas the supersymmetric gauge theory is 
formulated in only four space-time dimensions. Similarly, in developing the worldshcct description of A/" = 4 
supersymmetric gauge theories [17], we found it particularly convenient to add six extra dimensions, that 
is the index of took the values i = 1, 2, . . . , 8. The boundary conditions on the six new ^'s were strict 
Dirichlet conditions = on all boundaries, internal or external. At the same time wc added three new 
sets of b, c ghosts, which like the original set have strict Dirichlet boundary conditions. Since the extra g's 
and ghosts share identical boundary conditions, their contributions to the path integral exactly cancel: they 
are just dummy integration variables. 

For our purposes, to locally produce the necessary quartic counterterms in pure gauge theories, two extra 
dimensions and one extra set of b, c ghosts suffice. We thus have four transverse dual momenta corresponding 
to six dimensional space-time. Let us call the new dimensions r'^, k = 1,2 and wc can, if we wish, use a 
complex basis r^,r^. But here A,V do not represent helicity, but rather an analogous charge in the extra 
dimensions. Next we allow spurious with values ±1 of this charge to couple to two gluons as indicated in 
the top line Fig. 16. In order to guarantee that the spurion decouples, we insert a factor p^dr/da on the 




Figure 16: New cubic counterterms, with a spurion (dashed line) coupling to two gluons. 



spurion propagator near the interaction point. Because r*^ = on all worldsheet boundaries, the average of 
this factor over worldsheet fields vanishes, except when there is another such factor on the same time slice. 
In other words all the four gluon diagrams exchanging the spurion with internal propagators vanish, leaving 
only the contact contributions. By inspecting the coupling assignments shown in Fig. 16, we see that 

^s^r" = -{ad + bc)g' (147) 
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r^pV/ = -2(ac+6d)5^ (148) 
Altogether we have five adjustable parameters to produce two independent counterterms: 

r^AVV ^ ^AAVV^pAAVV ^ _(^ad+bc + Og' (149) 

r^vAv ^ ^AVAV^r^pV/ = -2(ac + 6d-C)/ (150) 

Since ^, a, b, c, d are arbitrary, there is more than enough flexibility to produce the necessary counterterms, 
and more generally to adjust them appropriately at each order in perturbation theory. It is perhaps most 
economical to employ the extra dimensions only to cancel the part of the anomaly due to UV artifacts, which 
would demand spin independence for this part: ad + be = 2{ac + bd). With this choice we then determine 
^ = -{ad-\-bc) = -1/1277^. 



11 Conclusion 

In this paper we have completed the lightcone gauge calculations of the scattering of glue by glue through 

one loop. Our results completely agree with those obtained using covariant methods [20]. In addition to 
obtaining the elastic amplitudes through one loop which are divergent in the infrared, we have also calculated 
their contributions to probabilities and have shown that infrared divergences cancel against contributions 
from extra gluons in the initial and final states. This is all in accord with the Lcc-Naucnbcrg theorem. 

The expressions for the final infrared finite probabilities, including the bremsstrahlung gluons, are ex- 
tremely compact, and are manifestly Lorentz invariant. Infrared divergences have been traded for a resolution 
parameter A characterizing unseen gluons and jets. 

All calculations were done in 4 space-time dimensions without the benefit of dimensional regularization, 
which means that (local) counterterms beyond wave function and coupling renormalization must be included. 
In spite of artificial p"*" = divergences (which raise the ugly possibility of requiring nonpolynomial coun- 
terterms), all necessary counterterms are polynomials in the external momenta of the degree dictated by 
power counting. We would like to underline here the fact that we did nothing sophisticated with p+ zero 
modes. We simply discrctizcd the integrals and omitted the zero modes [21,22]. Apart from coUincar 
divergences, which are only a problem for self energy insertions on external lines, this discretization pro- 
vides an apt infrared regulator for lightcone calculations. Our calculation shows that in all infrared safe 
calculations, including on-shcU gluon scattering with due care taken with jets and gluon bremsstrahlung, the 
continuum limit of this discretization is finite and yields all previously known results obtained in covariant 
gauges [18-20]. A sophisticated treatment of p'^ zero modes is not required. 

Finally we have discussed how to incorporate all of the counterterms we require in the lightcone worldshcet 
formalism. A particularly convenient way to do this is to interpret the QCD "string" dynamics as occurring 
in 6 dimensional space-time. We stress that the extra two dimensions are holographically generated on the 
"string" side of Field/String duality and are not present at all in the field theoretic description of the "field" 
side of the duality. Significantly, this can all be done while preserving a local worldsheet dynamics. 

For perturbative QCD, the next step is to prove that the lightcone gauge calculational procedure we 
have adopted carries through to all orders in perturbation theory. We believe it will because we have only 
needed to introduce strictly local counterterms, consistent with the concept that gauge violating artifacts are 
entirely associated with the ultraviolet part of the dynamics. The fact that we needed some counterterms 
that were not in the input classical Lagrangian is completely standard with the use of a gauge noninvariant 
regulation and should not obstruct the usual renormalization program. This is because the new counterterms 
introduced obey the power counting rules of renormalizable theories. 

If the renormalization program goes through as we expect, then the modifications we have made in the 
"bare" worldsheet description to accommodate the counterterms should suffice to all orders in perturbation 
theory. This would fulfill our ambitious goal of establishing a string theory dual of the large Nc limit of 
QCD, working entirely from the field theory side of the duality. Then the exciting prospect before us would 
be to use this duality to deepen our insight into nonperturbative aspects of the strong interactions. 
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A Bremsstrahlung Integrals 

In the evaluation of soft and coUinear bremsstrahlung cross sections we need to do several integrations over 
phase space. For the hard collinear case the integral was simple enough to treat in the text. Here we sketch 
the evaluation of the integrals for soft radiation which are more complex. 

At A^c = 00 each amplitude is the sum of gluon emissions from two neighboring lines, so the squared 
amplitude has two direct terms and a cross term. We first consider the transverse integration over the 
resolution of a direct term, which has the structure 



-^Direct(^^) = 



dk- 



(k 



- / kdk d(f) 
2 Jo Jo 



2kv cos < 



(A;2+t;2 + M2 + 2A;t;cos(;&)2 



(151) 



The (j) integral is easily done by transforming to a contour integral over z = e^'^ and evaluating residues. 
Changing variables from k to t = k"^ then leaves us with 



-^Direct (v) 



A'2 



dt 



[{t - t;2)2 + 2M2(t + t;2) + M4]3/2 



(A'2 



4 y'(A'2 - 1,2)2 + + 2A'2M2 + 2M2z;2 



4 



TT (M2 - t;2 + A'2 + _ y2y + + 2A^2j^^2 ^ 2M2t?) 

^2 2M2 



n A^ 

2 "^^^ 



V 



for A'2 > i;2; 



TT V 

2 t;2 _ A'2 



for A'2 < 



(152) 

(153) 

(154) 



where the last line applies as M — > 0. 

The transverse momentum integral of a cross term has the structure 



Ict 



-I 



dk 



{k + v)- (k + w) 



[(fc- 



/ kdk j 
Jo Jo 



2 + M2][(fc + t(;)2 + M2] 

+ vw cos a + kv cos ^ + kw cos(0 — a) 



d(j) 



(fc2 + t;2 + M2 + 2kv cos (j)) (fc2 + «;2 + M2 + 2kw cos(0 - a)) 



(155) 



Again the (p integral is easily done by converting to a contour integral in 2; = e^'^ . In this case we can 

simplify life by taking M = from the beginning, in which case the poles of the integrand arc at z = 
—vlk,—k/v,—we'^°'/k,—ke'^°'/w. The contour at the unit circle will enclose precisely two of these poles: 
z = —mm{v/k, k/v), — e'"min(w/A;, k/w), with four possibilities depending on the relative size of k, v, w. For 
definitencss let's assume that v > w. Then when w < k < v it turns out that the two pole contributions 
exactly cancel. Then the angular integral is 



1 



1 



vwe 
1 



+ 



vwe^' 
1 



fc2 — vwe fc2 — vwe^" 



ioTc k > V > w 

for V > w > k 
for V > k > w 



(156) 
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and we get with u = k'^^ assuming A'^ > v^^w^^ 

2 

1 1 



TT f 
"2/0 



du 



u — v'we~^°' u — vwe^' 



du 



1 



+ 



u — vwe~^°' u — vwe^' 



TT J A''* — 2A'^wu;cosQ! + w^w^ 

2 {v'^ + w"^ — 2vw cos aY 

TT , A'^ - 2A'2i; • w + v'^vP' 

= — In 71 1 

2 {v- w)^ 

IT ^ (v — w)^ 

--In^ ^ , 



for w"^ <v^ < A'2; 
for < A'2 < 



--In ^ , for A'^ <w^ <v^ 

2 V w 

Uv < w, the same result follows. The complete transverse integral for soft gluon radiation is the combination 

-^Total = -^DiroctC^') + -^Dircctl'w) — ^Cross 



— In ■ 



{V - W)^ 2 2 A/2 

A ,/i r, A ,9 „ „ , TT In — —5 , f or w < u < A ; 

2 A'4 - 2A'2t; • w + v^w^ M^e 



2 1,2 -A'2 ^2 M^e ' 

TT, (A'2-t;2)(A'2-t«2) 

= In- 



for < < 

for A'2 < w;2 < 



2 A'4 - 2A'2t; • w + t;2w2 ' 
For A;+ — > both u/A', w/A' — > so we see that the small region is insensitive to the resolution A'. 

B Evaluation of Bubble and Triangle Integrals 

We list here the q~ integrations of bubble and triangle integrands that are useful in analyzing coUinear 
divergences. First the six bubble integrands: 
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for < go+ < p+ (161) 
for < \q+\ < b+ I (162) 



— ^ 



—I 



—I 



dq 


1 




1 




for ^<qQ <V\ 


2-K 


mi 


2pf(9o- 


9cj"Pi/pf)^ ' 




dq^ 


1 




1 




for ^<qx <vt 


2-K 




2pJ(9i - 


9l"P2/P2)^' 




dq~ 


1 




-1 




for < < -vt 


2-K 


2 2 

9293 


2P3 (93 - 


93"P3/P3 )^ ' 




dq~ 


1 




-1 




for < g+ < 


27r 


2 2 

mi 


2^(90- 


itvjvtf ' 




dq~ 


1 






1 




2-K 


9o92 


2Pi2 [(9o 




)^ + 9o+(P^2-9o+)P?2/pi2'] 


dq~' 


1 






1 




2-K 


2 2 

9193 


2bi4l[(9i 
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Here we recall the notation qi = q — ki, where ki are the dual momenta, related to the gluon momenta by 
Pi = ki — ki-i, with ki = ko- We normally take k^ = 0. Next we list the triangle integrals: 



dq 

27r qoqtq2 



1 

2„2„2 
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^12 % for n+ < o+ < n+ 



2P^P^2(90 - Pl2 - - P]*2)P2/P2 )^[(90 - 9([Pl2/Pl2)^ + 9([ (P^2 " Qo)Pl2/Pl2} 

dq- 1 



27r q^qfqi 



for < g+ < p+ 



-it - Pi 

'^PtPlAiQo - qtPilPt?[{Q0+Pi - {dt +Pt)Pu/ptA? - {(lt+Pt){Qt - Pt)piA/ptA] ' 

for <qo < -pt 

dq- 1 

27r gigi^i 

Z^L 

2ptptAiQl - QtP2/ptmQl - dtPlAlptA? - (it + PtA)qtplA/pti] ' 

for 0<qt < -ptA 
'^^ +7^, for -Pu<Q+ <pt 



I 



^Ptpt ill - QlP2/ptf{Ql - P2 - {Qi - P2)P3/pt) 

dq- 1 



27r qoqhl 



for < (7+ < -p+ 



2^12^4 [(9o - 9+Pl2/Pl2)^ + 9o"(Pl2 - 90 )Pl2/Pl2^](gO " 9+^4/^4)^' 

go" - P12 

2^3^12(90 - P12 - {Qo - Pl2)P3^3 )^[(90 - g+Pl2/Pl2)^ + <lo{Pl2 " 9o")Pl2/Pl2^] ' 

for -pj <q^ < 



Once the g integrals have been performed, the transverse q integrals can be done using one or two Schwinger 
parameters to exponentiate the one or two denominators. 



C Evaluation of Box Integrals 

The box integrals we encounter can be most easily handled through the introduction of Schwinger parameters 
Ti,T2,T3, T4 for the internal lino propagators {q — fco)^^, — i/l ^ ^'2)^^, (<Z — ^3)"^ respectively. Since 

some of them are divergent in the ultra-violet, we also retain the worldsheet UV cutoff factors e"'^'' . The 
integration over q is then a Gaussian that is easily done by completing the square and shifting q ^ q + K , 
with 

TiA + S ' Ti4 ^ ^ 

where we use the shorthand T14 = Ti + T2 + T3 + T4. One then finds, using the Feynman parameters 
Xi = Ti/Tu that 

K16 -PiQ + Q'^Pi - X3K12 

K52 -Ptl + Q^P2 - XAK23 



X4K41 +pi — 

-LIA 
-114 



(164) 
(165) 
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SIC 

-K'ss p^q-q^p^+X2K23+xiKs4.-pt7^ (166) 

-tl4 

5K 

Kei pt Q - q~^P4 + ^3^ 34, + X2K 41 -ptTjr- (167) 

J 14 

and the Gaussian factor left over is just exp{—Ti4(x 1X3^^3 + X2X4Pi^) + 0{S)}. The 0{S) term in the 
exponent is negligible in the box integrals because the log divergence is insufficient to overwhelm it. Also 
when the Kij occur in the numerator of the box integrand the terms pfSK/Tu are negligible since they 
are 0(1) only when all Tj = 0{S) and they occur only in integrals convergent in this region. The Gaussian 
integration over q involves up to four powers of q as prefactors. 

/,^„^^^e-.--'-- = (168, 

Changing variables from the Tj to three of the Xi and Tx4 = T reduces the integration measure to d'^x5{l — 
Y^XijT^dT. In the first three cases it is safe to set (5 = 0, and the integral over T gives 



{xiX3pl2 + X2X4P\4Y ' 2{x\Xip\2 + X2X4Pi4) ' XiXapj^ + X2X4p\4 



(169) 



respectively. In the last two cases we must evaluate integrals that are log divergent for ^ ^ 0. We find, as 
(5^0, 

Jo WTW^' ~ -E^-7-ln(^^) 

f (f^e- ^ -^-.-ln(.^), f ^&^e-~-4-3,-31n(.^)(170) 

where 7 = — r'(l) is Euler's constant. 

We shall have use for the following combinations of momenta which arise in the box integrand after 
shifting q and sending ^ — > 

Ko = X2P1 + xs{pi +P2) - X4P4 (171) 
Ko-pi = X3P2 + X4{p2 + Ps) - xipi (172) 
K0-P1-P2 = X4P3 + xi (p3 + P4) - X2P2 (173) 

K0+P4 = XiP4 + X2{P1 + P4) - X3P3 (174) 

Finally, we list the integrals over Fcynman parameters that arise in the box diagrams. We use the shorthand 
notation (fix = Y^=i dxid{l — J2i Xi)- 

/■ ,3 w . ^^ 11 B\nB + A\nA AB / o , 2 ^\ 

L = jSxHxiX3A + X2X4B) = --+ ^^^^^^ +T2(:4TBFr+^'^ ^'''^ 



Li = J rf^x _ ^ " _ „ = ( TT^ + In^ ^ ) (176) 



1 _ 1 / 2 1jj2 ^ 

' X1X3A + X2X4B ~ 2{A + B) \ ^ B 



- \,X3A + X2X4B = 2(aTb) + 4iATW [ ^ b) ^'^^^ 



3 {X2,X4) ln{B/A) A / 2 , , 2 ^ 



= / \,X3A + X2X4B = 2{ATb) + A{A + Br V ^ b) ^'^^^ 
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D Evaluation of Eq(93) 

The numerator in the integrand, after introduction of Schwinger parameters and the appropriate shift of q 

can be replaced as 

K^,K^,K^,K^^ [p+g^ + X2X2I5 + a;ii^^4] [pt<l^ + ^^K^a + ^2K^i] 

X [pW + x:iK^2 + XiK^^Mq"" + XiK^^ + xiK^^] 
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+P^p|(9^<?'' - X2XapIJ2)[x2K!^:, + xiK^^][xaK'^^ + xiK^^] 
+Ptpt - xiX3plJ2)[x3K^^ + XiK^^][xiK^, + x^R-^^] 

+Ptpt{^''<l' - X2X4pIJ2)[xsK^2 + X4K^,][X3K^^ + X2KI] 

The second term in the first hne of the last equality came from the quantity 
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which can be greatly simplified using 
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Putting in the rest of the factors and doing the q integration yields the x integral 
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L2a{A,B) 



■'^23-'^41 _|_ -'^41-'^23 



Ptpt 



P2P3 



P1P2 



PiPt 



+ Lc{A,B) 
+ Lc{A,B) 



-^34-^41 + -^23-'^34 
ptpt 



+ 



-"■12^23 ^ -"^41 -'^12 



P1P2 



KX1KX4 



ptpt 



PiPt 



1 



167r2 



-(ln(5e^) + 2) - L{pI^,pIJ - 2p\^p\^Lab{p12.pIa) 



+Pu{ L2aa{Pu,p12) 



+P12' 



L2Aa{p12,p14) 



^23^^41 , -'^41-'^ 23 
+ + 
P1P2 

JV^2-'^.i4 



ptpt 



ptpt 



Pfpt 



+ Lca{p14,Pi2) 

+ ^Ca(p?2.Pi4) 



-'^34 -"-41 



ptpt 



K^,K^ ^ Kl2K^3 + ^41 ^^2 



-^23-^12 



-"■34^23 



ptpt 



+ 



PlPt 
-'^12-'^41 + -^^41-^ 34 

PxPt 
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The various L's in this formula arc Usted in Appendix B. Eq. 194 can be algebraically rearranged and cast 
into a nicer form. Putting in the coupling constant we find for the contribution of (93) 



if /(tree 
-g^2^AVAV 



1 



1 



2 2 
P12 



1 



2 3 
P12 



1 



2 4 
P12 



2 



H4 



16 (pf 4 + 2 



P?2 



1 



H2 



32b?4+P?2)' 



16(P?4+P?2)^ 



32(pf4+pf2)^ 
1 



P?2 



Pl2 



32 pU+pi^ 32ipU+pl^r 



1 ^12(^3" k2){k^ ^0^) _|_ 1 ^12(^3" k2){k, 



1 



48 (4 - fco+)p?4(fci+ - kt) 48 (4 - k+)pUkt - k+) 48 A;+ - k+ 

1 



48 ^2 



48/c^ - 48A;^ 



1 

+ 24 



fc+ - k+ 



2 /^'^ ^2 



Pl4('''3^ ^2^)^(^3^ ''^O^) Pl4{kt ^2^)^(^3~ ^0~) 



^14(^3' ^2~)(^2~ ^0~)(^l' ^2~) 
1 . o . „ J_1 

144 j 



^^2^^43) -^l0gP?4^e^ + 



E Complete Results of the Two Model Boxes 

There are basically two model integrands that has an adjacent pair of vertices with the same helicity. 



Assume for the moment that TZ = 1. These diagram can be reduced into triangles using Eq.(92). The first 
contribution is 



k^ <q+ < kt 

liq+-k+r(k+-k+){kt-k+) 



log{H,6e-y) + 1 (4 l^rikt kp{k^ k^) ^^^^^^^^^^ 



1(^ 
4 



(^3 ~ ki )pi4 



l{q+-kt){kt-q+)^^ 



4 {k+ - k+)p 



14 



k+ - 

klfif4 -k+Kf^^ -k+) 



/C3 Ki 



+ 



+ 



kt < ki 

1 (fcj" — fcj")(fc(}" — k^)^{k^ — A:^) ^ 1 {ki — — fc2^)(fct| ~ fci")(g"'" ~ ^1") 

1 (A:+ - kY){l4 - k^}{l4 - k- - kY + kY}(q+ - k+y 



{ko — k2){k'^ — k^') 



logiHsde-') 



1 {q+ kt){kt 9+) iog(/^,<5e-) + g+) j^Aj^v iog(/f,5e-) 



+ - 



4 (fc3+-fci+>?4 
l (fc3+-g+)^(fc+-fc2+)(feo+-fci+) 
8 /i^-j^ 
l(9+-fc+)2(fc3+-4)(fc3+-fc+) 



4 - fc+)p2 



log(ifdJe^) 



4 (^3+ - 



14 



^3 
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1^+ _ 4-+ 



+ - 



l (fc3+-g+)(g+-fc+) 
4 (fc+-fc+)p?4 

4 (fc3+-fc+)p?. 



l (g+-fc+)^(fc3+-fc+)(fc3+-fco+) 

8 /i/g fc"!^ 



Assuming its coefficient to be — C — A , ^ \ 
will be 



or -C-A 



(q+-k+r 



+ B 



(g+-fe+) 



log(if„,5eT) 



its contribution 



1 

8^' 



36 " 24 



log (p?2'5e^) 



(— fc;^ + kQ){—k^ + k'l){k2 k^){k2 — k^) 



+ 



+ 



w ^"^3 i^Q ^2 1 3 3 2 ^ 2 ^ \ 1 3 / i ^12 

24 24 



/C-j^ ~|~ fcg /Cq 

8" 



+ 74log4^T|^^^log4^ 
4Pi4 8 pii pU 



^14-^32 + infra red terms 



And it's impossible for this diagram to have poles at q'^ = kg or fc^. The reader can work out the infrared 
terms by partial fractioning the coefficient of each logarithm, and taking only the pole terms (but be very 
careful when applying these results to actual boxes, you might need to take a conjugation occasionally). 
The second case is very similar: 

-^n'^K'^ K'^ K'^ X 

g^2 5 -"■61-"-25-"-5,3-"-4,6 ^ 



+ - 



k^ <q+ < kt 

l (g+-fc2+)(fco+-g+) 
^ {ko ~ ^2^)^12 
l {q+-kt){k+-q+) 
4 {k+-kt)pl 



K-,KMHsSe^) + Q^kt kp{kt k+ , ^^^^^^^^^^ 



12 



+ 
+ 
+ 
+ 



fc^ < g+ < k'^ 
1 

4 



l{q+^kt){~q+ + k+). 



12 



1 {kQ ki)^{k^ k'^){k^ fcj") ^ 1 (^3" feo")(^3" ^2^)(feo^ ^]^)( g^+^o") 

4 



^3 ^'1 



K,3 fi,]^ 



1 (fc+ - A-+)(A:+ - A-+J(A:- - A-+ - k+ + A:+j(-r/+ + k+f 



4 {k^-kt)ph 
l {q+-k+){-q+ + k+) 

^ {ko ~ ^2^)^12 

l{q+-k+)i-q+ + k+) 



{k+ - k+)p 



l (g+-fc2+)^(feo+-fc+)(fc3+-fco+) 
8 k^ — fc^ 

l {-q+ + k+r{kt-kt)ik+-k}) 
8 k^ — k2 



\og{Hr5e') 



12 
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+ 



l {q+-k+){-q+ + k+) 
1 {q+ - k+){-q+ + k+) 



Assuming its coefficient to be -C-Aj^^^^+Bj^^:^ or -C-Aj^^^^+Bj^^^^, its contribution 



will be: 

1 

8^' 



36 ~ 24 



log {pl^Se^) 



{k+ - k+){kt - k+){k+ - k+m - fc„+) 



+ 



'-^ \ "'I "'3 T ft-i ft/Q T /S/Q "-3 T /ti ft'2 ^ "'2 "'3 •^"'0 "'2 / I ^ \"0 "'2 ^ i„„. I'l2 



24 



P?2 



24 



Pi4 



S + k^ + A:^ fc^ ^ 1 , ^12 I 1 fep" ^2^ , Pi2 



P?2 



H4, 



-fC^i i^Yo + infra red terms > 



And it's impossible for this diagram to have poles at q^ = k^ or fcg . 

F Complete Results of the Second Term of Eq.(104) 



1 (fcg" k2){k^ fc|j')((fcj~ k^)Pl4 ^ Pl2^3 Pl2^o) I 1 {^-3 ^2~)Pl4 + Pl2^0~ Pl2 



1 4, 



,2 Pl2 



log -^+7^ 
Pl4 



(4 - k+M.+pl^) 4 (fe3+ - kt)plM4+Pl2) 

1 Pl4Pl2{~^0 + ^3^)(^3~ ~ ^2')(^3~ — fco" + fc^ — fc^) 



^14-'^32 I 



16 ipl4+Pl2)Hkt-kt) 

1 P?2((-fc0+ + 24 - 4)P?4 +P?2(fc3+ - 4)) 



{kt-kt){pi+pi,rpi. 



^1A^32 



1 ^^12(^2 ^1 }(^^2 ^2 2A/^ /i/2 H~ ^2 "^^0 ^3 ^^3 ^1 ^1 ^0 ) 



16 



P?4 + P?2 



T— 7Z~i _ I , , r, ~. -fVqoiV, 



+- 



{kt-kt){pl,+Pl2) 

4 (4-4)(f?4+ff2)' 



(Pl4+Pl2)^ 



^32^A3 



1 

8^ 



1 k-^ 



1 



■^14-^^32 + Tc ('''s' ~ 4)^(~4 + '''2') 



8 p?4 ■ 16 



log 



2 P12 
Pl4 



16 



P14+P12 
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+ 



1 (^2 ^1 ) (^^2 ^2 ^2 2/C-|^ ^3 ^^3 ^1 ^^3 ^0 ^0 ^1 ) 



16 



± y nj-^ rvi^ \^ /vq fvt^ ''^Q A^g 3 3 1 3 2 1^ '^2 /-t"^ 



(4 



2 



'^32-'^ 43 



1 Pi2 kg + , l Pi2 (^3" k^ + ki ^2")^ r^A i>'V 

+ ~ 77jL , _i_s ^ o ^ TT:; 32-^^43 



16 



(P?4+f?2)' 



8 (fc3+-4)(p?4+p?2)3 



87r2 



1 Pi2(^J ki ){kt k^ -\- ki ^2")^ I 1(^3" A;^ + fci" ^2") 



16 



P14+P12 



+ 



8 (fc3+-fc2+)(p?4+P?2) 



-'^32-'»^43 



— (^3" k^)'^{k^ A;2~) + Yg^'^'S^ 'i'i~)(^3^ ^2~)(3fc3~ 2^2") 



+ g ^p2 ^■^14-^32 IQ^^S " kt)^^^^^ — — fc^) — - - 



k+ 

-"^14-'^32 



The results for other pole locatioins can be obtained from this by rotational symmetry and conjugation. 

If wc put these terms from each diagram together, there is usually some nice simplifications. 
For the fifth box of fig.7: 



■8^2^AVAV 



log^+TT^' 
Pl4 



1 



2 2 
P12 



2 3 
P12 



_'^{Pl,+Pl2? {Pl,+Pl2f 

And also a piece that doesn't fall into a tree: 



PI2 



+ log 2 



2 2 
P12 



'{Pl,+Pl2? 



PI2 



Pli+Pl2 



1 4 1 P12 

^9 log^x 



pIa 



2 Pi2( fep" + ^2") _|_ ^12(^3" ^2") _|_ 2 ^12(^3" ^2^) _ 2 ^^'^ kp^i ^0" + ^2") 



-^14(^0^ k't) Pui^^t ^t) ^14(^2" ^l") 



{k+-ktr 



^ ^ ( fcg + ) (2/03 fcg fc2 ) _ 2_ ^0 + '''2 



4-+ _ i'+ 

"'2 "'1 



+ 

+4- 
1 

8^ 
+ 



-4 



^14(^3" ~ ^2')(^3' ~ ko)iko ~ ^l) Pui^t ~ ^2')(^3^ ~ ^l")^ ^14(^3^ ~ ^2')(^3^ ~ ^o") 



-4 



2k^ k^ 



^2 



+ 4- 



(^3~ ~ ^0~)(^3^ ~ ^2~)^Pl4(^3~ ~ ^t) 



, (2^3~ ^2')Pl2 ^o" "I" ^2~)(2^3' ^o" ^2~) - kn + k^ 



ik+ - kt)p 



9^+ 1-+ 1-+ 
Zft,3 fi,Q fi,2 



1 



^14(^3" ~ ^2") (^2" ~ ^i") 
+ 2 



-"■32 -'^43 



«,3 fi,]^ 



+ 2 



^14(^3' kiYik^ fc^) ^14(^3' ^^)(^3' ^2^) 
For the fourth box of fig.7: 

5 



1 
8"^2 



" „2 /(tree 



Pi4 



P?2 



-'^ 32^43 



2 2 

P12 



2 3 
H2 



Pli 



PI2 



PU+P\2 '^{p\a+P\2Y 

2 



{pIa+p12? 
PI2 \ 



2 2 
P12 



PU+PU ipli+Pl2Y 

And also a piece that doesn't fall into a tree: 
1 



PU+PI2} 
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+ 



1 ^12(^3^ ^o")(^3" ^'2) _ ^Pui^'S ^■Q)i^''3 ^2) 



{^2 ~ ^o)p\A{^t K 



^2 



+ 



-(fc+-fc+)(fc+-fc„+)(fc+-fc+)p?2 

-214 + ^-(1 + ^-t 



Pui^t ~ ^o)i^^i ~ ^2) 

4 



kt 



^2 



^2 



(^2~ ^o')('^3' ^o')(^l' f^o)Pl2 



+ 4 



ik+-k+r{k+-k+)pi, 
1 



"■3 1^2 



-"-21-"- 14 



+ 4- 



9i-+ — h+ — 
"-S ~ n,2 



-2k+ 



^9. + ka 



{^2 kQ){k'^ '!'o~)Pl2 

1 



+ 



+8 



{^2 ^o')(^3~ 1^2) {kf 1^2 )p' 

1 

+ 8 



ko)ikt k2)Pi2 



J^21^i3 



-'^43-'^ 21 



k^ 



k+ 



^3" 



kt 



Puikt kQ){kt k^) Piiik^ ^2^)(^2^ ^o^)(^i^ ^2^) (k'^ ^o^)('^3^ ^2^)^14 



9i-+ — i-+ — i-+ 



{k^ — kQ){k^ fc2^)^Pi4 



-"■32-'^ 43 



i-+ — i'+ 



^(l)(^3" k2){kl k2)Pi2 



{kt-k+Y{kt-kt)v 



87r2 



5 X 



2 '^3 '^O o '^3 ^ 1 2 



■,Pl2(^3~ 



-'^43-'^32 



^o')(^3~ 



kt) 



kt ~ kt 



kr^ k> 



«-3 H.2 



+ 



+ 



+4 



1,+ 1,+ 
^3 ^2 



{kt kt){kt kt){kt kt)p^ 

- 4 



+ 4- 



^3 '^1 



Pi4iki — kt){kt — kt) 

{kt — kt){kt — Sfco" + 2fc3") 



12 (^2" kt)'^Pi2 
kt 



(kt 



kt)'{kt 



ko )Pl2 



^21-"^ 43 



kt 



ik'2 



kt)ikt kt){kt kt)Pi2 



+ 4 



{kt kt){2kt + kt 3kt) 



{kt kt)^Pi2 {kt kt)^{kt kt)p' 



12 



^43-*^ 21 



\ /i,2 j- ^3 ft-Q ^'1 '^'O ^1 ^3 ^^2 + ^0 ^3^0 ^^3 i tv^A tv^V 



kt){kt kt)^{kt ^o^)Pi4 



'^32-'^43 



For the seventh box of fig. 7: 



1 

8^2^ 



_^_„2 A tree 



And: 



1 .4i„„P?2 



87r2 
+4 



4 1 

5 log- 



PlA 



{kt 



kt){kt 



{kt kt)'^{kt kt){kt kt)p 



(ki 



kt)ikt kt){ kt 



2k1 



(2k 



kt){kt kt){kt kt) 



{kt kt)'^{kt kt)'^p\4ikt kt) 



-'^32-'^ 14 



{kt — kt)^{kt — kt){kt — kt)^Pi4 

^ jkt ~ kt){kt ~ kt) 

{kt — kt){kt — kt){kt — kt)'^Pi4 



14-'»-32 



^32-"-21 



-4 



{kt kt){kt kt) 



{kt-kt){kt-kt)^{kt-kt)p 



-A i^v I ^(^3" ^2')(^i' + ^3" '^kt){kt ^^l/^Aj^V 



^14^43+4- 



ftQ ) l^ft.]^ «,Q ) Pi2\l^2 '^0 J 



3 kt){kt kt){ 2kt + kt + kt) v , a 

-"■43-"-21 + ^ 



(fcg ^2 )(fcl ^2 ) 



{kt kt)'^{kt kt )'^Pi2ikt kt) 



{kt kt)^{k. 



'^21-^^43 



kt){kt kt)Pi2 



9^ X 



87r2 

I i~kt + ^3" ~ ^0" kt){kt — kt){kt — kt) 
pUkt - kt)-'{kt - kt)^{kt - fco+)2 
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Ti2 — h+'\3(l.+ _ l,+ \2(h.+ _ l.+ \2 

(/c2 A'Q ^0 ^0 ^1 ^^3 ^^3 ^1 ^2 ^1 ^2 )'^^32'^^14 

l'''3 "'0 ) ^^2 "■0 / V'*'! "-0 J P12 
( ^k<2 /^Q H~ fc]^ ^3 ^2 ^1 ^^3 ^0 ^3 ^1 )'^^21"^^43 

(fc.^ — fc(j")(fcj'" — fco")( — fc^ + ^3~ ~ + ^1") . 
fi,2 J I^Q ) {1^1 fi'2 ) P12 

(3^2 ^2 ^^3 ^2 ^^2 ^0 ^0 ^1 ^3 ^1 ^3 ^0 )"^^43'^'^21 f 



G Details on Triangle-like Diagrams with Collinear Divergences 

After ascertaining that the combination of 

A \{-kt + k+){-k+ + q+) , {-k+ + k+){k+-q+) 



2 (<z+ - ktr 



2 2 9 ~^ 2 2 2 



doesn't have coUinear divergence, we simply hst its contribution here. Note in the results below, we have 
put back in the term 

A Kl,K^^^ (^kt + k+){q+ - k+) + {k+ - k+){q+ - k+) Kl^K^ K^^Kl^ 

2 (9+ - kt? qlqhl {q+ - k+f qlq^hl 

(more specifically, the above contains all the triangle-like terms in Eq.(104) and the two subtractions from 
the first term of Eq.(104)). 

1 4 
A;^ < 9+ < kt 

8 (q^ - ktnk^ - k^) '"^^""^^ + 4 (,+ - kt){kt - kt)pUkt - k}) '''''''' '^^^""^^ 
I (fc3+-g+)(g+-fc2+)(fc+-fc+) 1 fc3+-g+ K^K^^ndUA 



fc+ < g+ < k+ 

1 (g+ - fc2+)(-g+ + k+){k+ - fc+)(fc3+ - k+) 1 1 

lOg(Jis) + -7^-^ x 2 -^43-^32 iOg(iisj 



(A;o+-fc2+)(g+-4)2 ''^ 4{k+-k+)p 



14 
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+ - 



1 (kt - kt){kt - q+){k+ - k+) 

{q+-kt){k+ ^kt) 
l (g+-fc+)(fc3+-g+)(fco+-fc+) 
8 {q+ - fc+)2 

1 ('1^ -k^)ij4 -^^t) 



log(Fd) + 



k+-q- 



4(fc+-fc+)p2^(g+-fc+)(fc3+ 



1 (fc+-fc+)(fc3^-g+)(fco^-fc+) 1 



8 (g+_fc+)(fc3+-fc+) 

k^ < q+ < k2 



4 {k+ - k+){q+ - kt)p 



8 iq--ktnk^-kt) - ~^pUk^-kt)iq--kt)''^'''^' '''^'''^ 



4(fco+-fc+K4 



+ 



fc+ 



1 iq+-k+)ik+ 



4(g+-fc+)(fc3^ 



+ - 



^{k+-k+){q+-kt)pl/^^^ 



log(i?d) 

-k+) 



4^4^"' " '^ '^ 4;,f,(g+ - kt){kt - ktm - kt) 



Again, we can integrate out whatever can be integrated out, and sweep the rest underneath the infra red 
term rug. Thus we have: 

k3){kQ — k^){k^ — kQ) ^ 1 ki — k^ 



1 4 , Pl2 



1(4 



The combination of: 

A j^5^^3.2Pl4 ^ ^ ^^5^^4,3^14 

2 qhl<ipl2 



2 9o9l92Pl4 



+ 



(-fco+ + fc3+)(-fco++g+) , (-fc+ + fc+)(fco+-g+) 



2 2 2 

9o9293 



2 2 2 



2 (<z+ - k+Y 

A ik+ - k+){q+ - kt) + {k+ - k+){k+ - q+) K^,2Kl, Kl,K^^, 

^ ' ihhl {q^-kt?qlqlqlql 



2 (9+ - kt) 

This is the rotation by 90 degrees clockwise of the previous one. Its contribution is: 



^ 4 
8^^ " 

k-^ <q+ < kt 

1 (_fc+ + + _ ^^^^^^^ 



1 

+ - 



(-A=o+ + 9+)(-fc2++4) 

-kt + q'^ r^V 



— if3^2Xl'3log(/7,) 



^i.-kt + q+)pl2{-kt + kt) 
1 (-fc+ + g+)(-fc+ + fco+) 



i^^2i^4^3l0g(if,) + - 

i^3^2^4Xjlog(iI„) 



^{^kt + q+)pU-kt + kt) 
1 + fco+)(-fc+ + k+){-kt + g+) 



log(ff.) 



4 (-fc+ + - kt)pi,{-k+ + fc+) 

1 {kt - kt){-kt + g+)(-fc2+ + kt) .(jj ^ 

8 ^k+-kt){-kt + q+) 
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, 1 (-fc+ + fco+)(-fc+ + g+) .^A..Vi 

fc;^ < g+ < A;^ 

1 (g+ - fc+)(fc+ - k+){-k+ + fc+) 1 (-fc3+ + kt){-kt + q+)(-kt + q+){kt - k+) 

8 (_fc+ + 9+)(-4 + /c+) 8 (_fc3+ + fc+)(-A;+ + g+)2 *^ 

1 ..A .(rr s l (g+-fe+)(fc+-fc+)(-fc3+ + g+) 

" 4(-4+fc^).f.(-fco^+.-^)''^^''-^°^^'''^^-8 i-k^+,-r '''^'''^ 

1 (-fc+ + fc+)(-fc+ + g+) 1 (fc+ - fc3+)(-fc+ + g+)(-fc+ + fc+) 

^ 4 (-^3+ + kt)pUkt - fc+)(-fc+ + ) ^32^43 l0g(^^.) 



1 (-4 + kt){-k+ + k+) ^ 1 {k+ - k+){-kt + k+){-k+ + q+) 



4 -k+ + k+ 8 (-k+ + q+){-k+ + k+) 

1 fco'fcj'" — 2k^ki — 2fcJ"(7+ + g+fcj^ + k'^k^ + l^^i j^a^j^v^ log{H ) 



log(i^«) 



12 



4 (A;+-4)(-fco++g+)(fco+-fc3+)F 
1 {-k+ + k+){-k+ + q+) yA.^Vi 1 -4 + g+ K^K^^nr(ff\ 

1 (-fc+ + g+)(fc+~fc+)(-fc3++g+) -kt + q+ 

S ^ 7+ Ty2 iOg(//;j + - +^ 2 / ;+ , ,+n -^32-^43 iOg(-H;) 

1 (_4 + fco+)(-fc3+ + + q+) 



i-k^ + q+){kt - k+) 



+ 4 - fc+K2(-4 + <I^){-kt + k}) '''''''' '^^^"""^ 
Its finite contribution is: 

87r2^ '°S?4 



/Cq^ k2 4 (/Cq" )p^2 



The combination of: 



gives: 



2 g|9o9?P?2 



? 9 ? ? 



k+r 



{-k+ + k+){-k++q+) ^ (-fc+ + fc+)(A-+-g+) 



9 9 9 

mm 



9 9 9 

mm 



2 (g4 

A i^r,6^5^ (fe3''-fe2+)(g+-fco+) + (feo^-fe3'")(4-9+) , K^K^a K^,5Kl^ 

2 (9+ - k+r ihlQl {q+ - k+r qhhlql 
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+ 
+ 



1 + q+)(k+ - k+){-k+ + k+){-k+ + q+) 



{-k+ + k+){-kt+q+Y 
2^ + *'' K^^K)[^\og{H,) 



log(F«) 



1 ~kt + g+ 

4(-fc3+ + (Z+)(-fc3+ + fc2+)Pl4 



2 



+ i:^jz^r^f^2K,,iog{H,) + - -^^^ iog{H,) 

1 {-kt+k+){q+-k+) 1 (-fc+ + q+){-kt + fc+)(9+ - 4) i„ .ff ^ 

~^pU-kt + kt){-kt + g+) '''''''' ip?4(-4 + <Z+)(fc0+ - + ^32^^43 l0g(i^.) 

< g"^ < k^ 

1 (-fc3+ + kt)iq+ - kt){-k+ + g+)(fc+ - fc+) 1 fe+-fco+ 

" 8 ^ ^(-fc^+A;^(-A^+g+)^ + i(-fc^+4)P?4(4-^^)'^^'^^ ^^'^""^^ 

^ 8 (-fc3^ + fc^)(-fc|+g+) + 4(-fc+ + fc,-)p?,(fc^-4)^32^- 

l (_fc+ + fc+)(_fc+ + g+) 1 {-k+ + q+)i~k+ + k+) 
~ 8 ;+ I 4- ^Og[H4) - - ,.2/7+ T+Tr-f^—n^^ 32^ 43 iOg(-Hdj 

l (-fc+ + fco+)(-fc++g+) 1 -fc+ + g+ r^AyV, 



(-fc3+ + g+)2 4pl^i^kt + kt){-kt + q+y 



{-k+ + q+){-kt+kt) ^A^^v 



4 (-^3+ + q+)pUk+ - k+K-kt + k+) 

k^ < q~^ < k^ 

1 (-fc3+ + fc + )(-fc+ + k+K-k+ + g+)^ 
8 (-fc3++9+)^(-4 + fco+) 

^ 4 - fc3+>f,(-fe3+ + g+)(-4 + 4)^-i^43l0g(i^.) 

l (_fc+ + fc+)(-fc+ + g+) 1 (-fc+ + g+)(-fc3+ + fc+) 

~ S ^+ I + ^OS{Hd) - T77+ ■+N 2 / L+ , +N/ ,+ , t.+ A -^32-^43 l0g(i?d) 

8 -kj+q+ ^[k^ -ksjptii-k^ + q+){-kJ + kj) 

l (-fc+ + fco+)(-fc++g+) 1 -fc+ + g+ r^A^^V, 

, 1 (-fc^ + g+)(-fc3+ + fc+)(-fco+ + g+) 1 -fc^ + g+ 



Its finite contribution is: 



1 4 , H2 

—5^ log 

8772 



1 (fcj" fc3")(fc]^ k2){k^ fco") ^ 1 fcj' fc.^ 
8 kt-kt ^{kt-k+)pl ' 
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The combination of: 

. A K^Kl^ \ {-k+ + k+){-k+ +q+) {-kt + k+m-q+y 

A K^^.Kl, (fc+ + fc+)(g+ - kf) + {kt - kt){q+ - kt) K^Kl, Kl,K^^, 
+ 2 (5+ - k+r q\qlql {q+ - k+f 



2;ives: 



+ 



kt <q^ < kt 



^ y ^ "-3 

log(ff.) 



1 {kt - kt){k+ - kt) ^ 1 {kt - kt){kt - kt){kt - q+) 



kt-kt 8 {kt - q+){kt - k+) 



, 1 1 ji'A Tv^v rr \ 1 kt kt 2kt kt ^q~^kt + q~^kt + q~^kt + ktkt 

+ 4(fc+-4K2 " " " 4 {kt^ktm-q^){kt-kt)v\, ''^^''^^ 

1 (fco+ - fc+)(fc+ - fc3+)(-g+ + kt) 

~ ~Hkt-q^){kt^kt)pl/'''''^ '''^'''^ - 8 {kt~kt)ikt-q-) '''^"'^ 

1 ikt-kt){-q++kt) r.Ar.V.(rr. 1 {kt - kt){kt - kt){kt - q+) 

^ ~Hkt-q-)pUkt-kt)^'''''''' - 8 {kt-q^){kt-kt) 

1 fco+ - g+ J., , . _ Hkt-kt){kt-<l^){kt-Q^) ,.^ . 

4(fc+-fc+)(fc+-g+)pf2 

kt < Q'^ < kt 

1 (fco+ - q+){kt kt){kt kt) 1 (^3+ - - + kt){kt - kt) . 

8 {kt - q+){kt - kt) ^^S^^^'^ + S (fc+_fc3+)(fc+_,+)2 

^if,3^21 log(/^0 - 8 {kt-q+){kt-kt) 

^43^21 l0g(i?«) - g (fc+ - log(^r) 



^pUkt-q+){kt 


~kt) 


1 /,:+ - q- 




^ {kt ~ kt){kt — 




1 kt - q+ 





2 



kt < q^ < k. 

8 ^ (fco^-fcJ)(fc^-,+ ) ^"^^""^^ - -Hkt -q^){kt -kt)pl/'''''^'°'^'''^ 

1 - kt){kt — kt){kt - q'^) , ^ ^ 1 (fc(j' — fcj')(fc^ — g+) ^ ^ i„„/-zj \ 

8 ) " i(4-.+ )(4-fc3;^)M2 " 

IT^i^^^TT^^T+V^^^^- + ^ (^0^ 'S^^''{-^P^%7 ''^^ iog(g.) 

4(A;J -g+)(A;J -fcj)pf2 8 {kt - q+){kt - kt) 
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1 k+ - q+ 



'12 



Its finite contribution is: 




Pi 2 1 (fco" k^){k^ fc(^)(fcj'" fc2") 1 k^ ^2 



'-21 



[8 fco+ - k+ 4 (A:o+ - 



'12 
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